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This  thesis  develops  new.  modified  Anderson-Darling  and  Crainer-von  Mises  goodness- 
of-fit  tests  using  the  mean  rank  and  median  rank  plotting  positions  for  the  normal  distri¬ 
bution  with  specified  parameters.  The  complete  critical  value  tables  are  presented  for  each 
test.  These  tables  can  be  used  to  t*>st  whether  a  set  of  observed  values  follows  the  normal 
distribution.  Additionally,  the  power  tables  are  presented  for  each  distribution  to  help  the 
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Abstract 

This  thesis  improves  the  powers  of  the  Anderson-Darling  and  Cramer-von  Mises 
goodness-of-fit  tests  using  the  mean  rank  and  the  median  rank  plotting  positions.  A 
Monte  Carlo  simulation  of  5000  repetition?  is  used  to  generate  the  critical  values  from 
n  —  4  through  n  =  80  for  the  significance  levels  of  a  =  0.20,  a  =  0.15,  a  =  0.10,  a  =  0.05, 
and  a  —  0.01.  An  extensive  power  study  is  performed  to  compare  the  powers  of  the 
modified  tests  to  the  known  tests  for  the  uniform,  the  exponential,  the  double  exponential, 
the  lognormal,  and  the  Weibull  distributions. 

The  powers  of  the  A-D  tests  modified  by  the  median  rank  plotting  position  and  the 
mean  rank  plotting  position  are  approximately  0.01%  higher  than  the  unmodified  tests. 
The  power  of  the  C-VM  test  is  10%  higher  than  the  unmodified  C-VM  test  for  the  Weibull, 
the  lognormal,  the  exponential,  and  the  dr  >le  exponential  distributions  when  the  median 
rank  plotting  position  is  used.  Using  th  an  rank  plotting  position  also  improves  the 
power  of  the  C-VM  test  for  the  uniform  ar  Weibull  distributions;  while  it  reduces  the 
power  for  the  exponential,  the  double  expont.^., .  ,  and  the  lognormal  distributions. 

Finally,  both  A-D  and  C-VM  tests  modified  by  using  the  median  rank  plotting  posi¬ 
tion  have  more  power  than  the  tests  modified  by  using  the  mean  rank  plotting  position. 
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NEW  MODIFIED  ANDERSON-DARLING  AND 
CRAMER- VON  MISES  GOODNESS-OF-FIT  TESTS 
FOR  A  NORMAL  DISTRIBUTION 
WITH  SPECIFIED  PARAMETERS 


I.  Introduction 


1.1  Background 

The  problems  facing  the  Air  Force  continue  to  grow  in  size  and  complexity.  Simu¬ 
lation  modeling  is  one  of  the  most  commonly  used  techniques  to  resolve  these  problems 
(16:2).  Decision  makers  evaluate  the  proposed  solutions  by  using  the  results  of  the  simu¬ 
lations,  and  in  turn,  are  more  able  to  plan  future  activities  more  efficiently. 

Collection  of  data  is  one  of  the  most  important  phases  of  the  simulation  process 
(16:43).  After  the  collection  of  the  data,  analysts  compare  this  data  to  a  theoretical 
probability  distribution.  Then,  the  problem  facing  the  analysts  becomes  how  well  the  data 
fits  the  hypothesized  distribution.  The  test  to  check  if  the  hypothesized  distribution  fits 
the  data  is  called  a  “goodness-of-fit”  test.  A  goodness-of-fit  test  measures  the  degree  of 
agreement  between  the  distribution  of  an  observed  data  sample  and  a  theoretical  statistical 
distribution  (5:187).  If  such  tests  give  a  good  fit,  then  the  hypothesized  distribution  can 
bo  used  in  simulation  modeling. 

The  most  commonly  used  goodness-of-fit  techniques  are  the  Chi-square  and  the 
Kolmogorov-Smirnov  (K-S)  tests  (5:168).  Other  goodness-of-fit  techniques  include  Anderson- 
Darling  (A-D)  and  Cramer- von  Mises  (C-VM)  tests  (6:109-114).  These  tests  have  been 
modified  several  times  in  order  to  get  a  better  goodness-of-fit  for  the  normal,  uniform. 
Laplace,  exponential,  and  Cauchy  distributions  using  known  and  unknown  parameters. 
The  latest  modification  for  the  A-D  and  C-VM  goodness-of-fit  statistics  was  performed  by 
Woodruff  for  the  logistic  distribution.  He  derived  the  critical  value  tables  for  the  logistic 
distribution  with  unknown  shape  and  location  parameters.  His  power  study  indicated  quite 


1-1 


good  power  against  uniform  and  exponential  alternatives.  The  modified  A-D  and  C-VM 
tests  had  more  power  than  the  modified  Kolmogorov-Smirnov  test  (19).  However,  not  all 
distributions  have  been  successfully  examined  for  goodness-  of-fit  when  the  parameters  of 
the  distributions  are  unknown. 

The  normal  or  Gaussian  distribution  is  the  most  important  distribution  in  probability 
and  statistics  (16:32).  Much  of  the  data  collected  can  be  described  by  a  normal  distribution 
or  can  be  said  to  be  normal  as  a  first  approximation.  After  the  data  is  collected,  one  can 
find  that  the  data  is  often  approximately  normally  distributed.  Further,  the  assumption 
of  the  normal  distribution  as  the  basis  for  analyzing  the  sample  data  is  often  used  in  the 
development  of  advanced  statistical  theory.  It  can  be  concluded  that  the  distributions  of 
many  sample  statistics  tend  toward  normality  as  the  sample  size  increases. 

The  investigation  for  a  better  goodness-of-fit  test  which  can  be  applied  to  the  normal 
distribution  still  continues.  No  tests  have  been  developed  for  the  normal  distribution  with 
known  parameters  by  using  the  new,  modified  Anderson-Darling  and  Cramer-von  Mises 
statistics  to  determine  a  better  goodness-of-fit.  This  test  would  be  useful  for  the  Air  Force 
so  that  analysts  can  determine  whether  a  random  sample  of  data  taken  from  an  observed 
phenomenon  behaves  according  to  the  normal  distribution. 

1.2  Problem  Statement 

Bow  can  the  existing  Anderson-Darling  (A-D)  and  Cramer-von  Mises  (C-VM)  tests 
be  modified  to  produce  a  new,  better  goodness-of-fit  test  which  can  be  applied  to  the 
normal  distribution? 

1.3  Research  Objirtives 

The  objectives  of  this  thesis  are  to: 

L.  Generate  and  document  the  known  A-D  and  C-VM  critical  value  tables  for  the  normal 

distribution. 


2.  Generate  and  document  the  A-D  and  C-VM  critical  value  tables  for  the  normal 
distribution  with  the  modified  A-D  and  C-VM  statistical  values.  These  tables  will 
be  used  to  test  goodness-of-fit  when  the  parameters  of  the  distribution  are  known. 

■3.,  Conduct  a  power  comparison  between  the  modified  and  unmodified  goodness-of-fit 
tests  to  determine  which  test  can  best  detect  a  false  normal  distribution  hypothesis. 

1.4  Methodology 

The  research  effort  will  consist  of  the  following  steps.  The  first  step  is  to  produce 
A-D  and  C-VM  critical  value  tables  for  the  normal  distribution.  The  second  step  of  the 
research  is  to  produce  the  modified  A-D  and  C-VM  critical  value  tables  for  the  normal 
distribution  with  the  new  modified  A-D  and  C-VM  statistical  values.  The  final  step  will 
be  to  compare  the  powers  of  these  two  modified  tests.  For  each  step,  computer  subroutines 
will  be  written  as  required  to  perform  any  necessary  calculations. 

1.4.1  Critical  Value  Tables.  During  the  first  step,  critical  value  tables  will  be  gener¬ 
ated  with  the  known  A-D  and  C-VM  statistical  values.  Random  deviates  for  a  fixed  sample 
of  size  n  will  be  generated  from  a  normal  distribution  using  the  Monte  Carlo  simulation. 
To  evaluate  the  effect  of  sample  sizes  on  the  critical  values,  sample  sizes  of  4  through  80 
will  be  used. 

The  estimated  parameters  will  be  used  with  the  ordered  normal  deviates  to  calculate 
the  values  of  the  hypothesized  distribution  function  at  the  order  statistics.  Based  on 
the  hypothesized  and  sample  distributions,  the  A-D  arid  C-VM  statistics  will  be  calculated 
using  the  appropriate  equations.  Each  of  these  steps  will  be  repeated  5000  times  to  generate 
5000  independent  A-D  and  C-VM  statistical  values.  These  critical  values  will  be  used  to 
obtain  the  critical  values  of  the  tests. 

The  same  procedure  described  in  step  one  will  be  repeated  in  step  two  using  the 
modified  A-D  and  C-VM  statistical  values.  Then,  two  complete  sets  of  critical  values  will 
be  presented  for  the  modifu  1  and  unmodified  goodness-of-fit  tests. 
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1-4-2  Power  Comparison.  Then  the  research  will  compare  the  powers  of  the  mod¬ 
ified  A-D  and  C-VM  tests  to  determine  which  tests  can  best  detect  a  false  normal  distri¬ 
bution  hypothesis.  The  power  of  a  statistical  test  is  the  probability  of  correctly  rejecting 
a  false  hypothesis.  The  null  hypothesis  that  a  set  of  sample  deviates  follows  a  normal  dis¬ 
tribution  will  be  tested  against  the  alternative  hypothesis  that  the  sampm  deviates  follow 
some  other  distribution. 

The  A-D  and  C-VM  test  statistics  will  then  be  calculated  under  the  null  hypothesis 
that  the  random  deviates  follow  the  normal  distribution  with  the  specified  parameters. 
The  calculated  statistics  for  each  test  will  then  be  compared  to  the  corresponding  critical 
value  tables  obtained  in  step  one  and  step  two  to  determine  whether  to  reject  the  null 
hypothesis.  After  repeating  this  oOOO  times,  the  number  of  times  each  statistic  exceeds 
the  respective  critical  value  will  be  counted  for  each  sample  size.  This  total  will  be  the 
number  of  times  the  null  hypothesis  has  teen  rejected.  The  power  of  the  test  for  each 
alternative  distribution  will  be  the  number  of  times  the  null  hypothesis  is  rejected  divided 
by  5000.  This  power  comparison  will  be  performed  with  some  alternative  distributions. 
These  alternative  distributions  are  the  uniform,  the  exponential,  the  double  exponential, 
the  Weibull,  and  the  lognormal  distribution... 


II.  literature  Review 


2. 1  Introduction 

A  literature  search  has  been  completed  on  the  goo-  ■  f-fit  techniques.  Research 
studies  published  and  unpublished  have  been  investor  "•  .  -r.^ugh  the  Current  Index  to 
Statistics,  the  Defense  Technical  Information  Center  t  -  i r  a..J  previous  AFIT  theses 
efforts.  The  following  paragraphs  will  review  the  literat  -  pertinent  to  this  research 
proposal.  Specifically,  the  discussion  covers  the  topics  of  h  .  '  .  .jsis  testing,  test  statistics, 
and  power  studies. 

2.2  Discussion 

From  the  earliest  days  of  statistics,  statisticians  have  begun  their  analysis  by  esti¬ 
mating  a  distribution  for  their  observations,  and  then  they  have  checked  on  whether  or 
not  this  distribution  is  true.  For  years,  a  vast  number  of  test  procedures  have  appeared, 
and  these  procedures  have  come  to  be  known  as  goodness-of-fit  tests  (6:1). 

The  goodness-of-fit  test  determines  whether  or  not  sample  data  could  have  been 
generated  according  to  a  particular  distribution  function.  It  achieves  this  by  comparing  the 
sample  frequency  distribution  to  the  assumed  theoretical  frequency  distribution  (12:597). 

Occasionally,  assumptions  are  made  regarding  a  distribution  when  there  rp&y  be  little 
prior  statistical  information.  In  these  cases,  it  is  more  desirable  that  there  be  verification 
that  the  distribution  fits  the  data.  Comparisons  of  test  data  with  an  assumed  distribution 
are  accomplished  by  a  goodness-of-fit  test.  Such  a  test  evaluates  how  well  the  data  fit  the 
assumed  distribution  or  determines  if  there  is  evidence  of  disagreement.  The  goodness-of- 
fit  tost  is  for  disagreement  rather  than  for  agreement  of  the  distribution  with  the  data. 
The  data  are  assumed  to  fit  the  distribution  unless  there  is  sufficient  evidence  to  disprove 
the  assumption  (1.-72). 

One  potential  disadvantage  of  goodness-of-fit  tests  is  that,  even  though  an  assumed 
fits  the  data,  there  may  be  other  distributions  which  fit  the  data  equally  well 
or  perhaps  even  better.  Another  distribution  may  be  of  a  different  type  or  it  may  be 


the  same  type  of  distribution  but  with  different  parameter  values.  Consequently,  when 
applying  a  goodness-of-fit  test,  it  is  essential  to  remember  that  the  test  only  rejects  the 
assumed  distribution  when  there  is  definite  evidence  that  the  distribution  is  incorrect;  it 
does  not  mean  that  the  selected  distribution  is  the  best  one  when  there  is  not  sufficient 
evidence  to  reject  it  (1:73). 

2.2.1  Hypothesis  Testing  and  Test  Statistics.  Hypothesis  testing  is  the  process  of 
inferring  from  a  sample  whether  to  accept  a  certain  statement  about  the  population  (5:76). 
Any  statistical  test  of  hypothesis  works  in  exactly  the  same  way  and  contains  the  following 
steps: 

1.  The  hypotheses  are  stated  in  terms  of  the  population. 

2.  A  test  statistic  is  chosen. 

3.  A  rule  is  made,  in  terms  of  possible  values  of  the  test  statistic,  for  deciding  whether 
or  not  to  accept  or  reject  the  null  hypothesis. 

4.,  On  the  basis  of  a  random  sample  from  the  population,  the  test  statistic  is  evaluated, 
and  a  decision  is  made  to  accept  or  reject  the  null  hypothesis  (14:429). 

The  hypothesis  to  be  tested  is  called  the  "null  hypothesis.”  and  denoted  by  The 
alternative  hypothesis,  denoted  by  Ha.  is  the  negation  of  the  null  hypothesis  (5:77). 

The  functioning  parts  of  a  statistical  test  are  the  test  statistic  and  an  associated 
rejection  region  '/>:78).  The  test  statistic  is  a  function  of  the  sample  measurements  upon 
which  the  statistical  decision  will  be  based.  If  for  a  particular  sample  the  computed  value 
of  the  test  statistic  falls  in  the  rejection  region,  the  null  hypothesis  H0  is  rejected  and  the 
alternative  hypothesis  Ha  is  accepted.  If  the  value  of  the  test  statistic  does  not  fall  into 
the  rejection  region,  then  H0  is  accepted  (14:429). 

There  are  two  ways  of  making  an  incorrect  decision  in  hypothesis  testing.  These  are 
cit  her  rejecting  a  true  null  hypothesis  or  accepting  a  false  null  hypothesis.  These  two  error 
types,  type  I  and  II  respectively,  have  associated  with  them  certain  probabilities  of  the 
errors  being  made  (5:79). 
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2.2. 1.1  Chi-square  Test  Statistics.  The  oldest  and  best,  known  goodness-of- 
fit  test  is  the  Chi-square  test.  Karl  Pearson  abandoned  the  assumption  that  biological 
populations  are  normally  distributed,  introducing  the  Pearson  system  of  distributions  to 
provide  other  models.  The  need  to  test  fit  arose  naturally,  and  in  1900  Pearson  invented 
his  Chi-square  test.  Modern  developments  have  increased  the  flexibility  of  Chi-square 
tests,  especially  when  unknown  parameters  must  be  estimated  in  the  hypothesized  family 
(6:63).  The  Chi-square  test  compares  observed  frequencies  with  expected  frequencies  of 
the  hypothesized  distribution  function.  It  is  restricted  to  large  samples,  approximately  25 
or  greater  (14:310).  The  major  advantages  of  this  test  are  that  it  can  be  applied  to  discrete 
populations,  and  it  can  be  modified  when  parameter  values  are  unknown  by  reducing  the 
number  of  degrees  of  freedom;  whereas  tiie  Kolmogorov-Smirnov  test  cannot  be  applied  in 
these  situations  (13:68). 


2.2. 1.2  Statistics  Based  on  the  Empirical  Distribution  Function  (EDF).  A 
general  class  of  statistics  used  for  the  goodness-of-fit  test  is  called  empirical  distribution 
function  statistics.  The  empirical  distribution  function  is  a  step  function,  calculated  from 
the  sample,  which  estimates  the  probability  distribution  function.  EDF  statistics  are 
measures  of  the  discrepancy  between  the  EDF  and  a  given  distribution  function  (6:97). 
For  a  given  random  sample  of  size  n,  Let  be  the  ordered  statistics;  suppose 

further  that  the  dismbution  of  X  is  F(x).  Then,  the  empirical  distribution  function  is 
/•’„{. e)  defined  by 

_  ,  ,  Number  of  observations  <  x  ^  . 

Fn(x)  = - —  (2.1) 

n 

where 


-  oc  <  x  <  oo 


More  precisely,  the  definition  becomes  (6:98-99) 


Fn(x)  =  0, 

*  <  XH) 

(2.2) 
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Fn(x)  =  1. 
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Thus  Fn(x)  is  a  step  function,  calculated  from  the  data;  as  x  increases  it  takes  a  step 
up  of  height  £  as  each  sample  observation  is  reached. 

The  Kolmogorov-Smirnov  statistic  (K-S),  which  is  one  of  the  most  well-known  EDP 
statistics,  is  the  largest  vertical  distance  between  the  completely  specified  hypothesized 
cumulative  distribution  function  (CDF)  and  the  observed  empirical  distribution  function 
(6:204).  In  some  cases,  the  analyst  may  wish  to  specify  completely  an  entire  distribution 
and  then  test  whether  or  not  his  sampled  population  has  exactly  that  distribution  in  every 
detail.  The  K-S  tests  are  appropriate  in  this  situation,  if  the  entire  specification  can  be 
made  prior  to  sampling  (2:296). 

There  is  a  controversy  between  the  K-S  test  and  Chi-square  test  about  which  test 
is  more  powerful,  but  the  general  feeling  seems  to  be  that  the  K-S  test  is  probably  more 
powerful  than  the  Chi-square  test  in  most  situations  (2).  D5Agostino  and  Stephens  showed 
that  EDF  statistics  were  usually  much  more  powerful  than  the  Pearson  Chi-square  statistic 
(6:110). 

Two  other  goodness-of-fit  tests  based  on  EDF  statistics  are  Anderson-Darling  (A-D) 
and  Crainer-von  Mises  (C-VM)  tests.  The  idea  for  the  C-VM  test  is  based  on  the  squared 
integral  of  the  difference  between  the  EDF  and  the  distribution  being  tested.  The  C-VM 
test  statistic  is  defined  as  (6:101): 


C  -  VM  = 


1 

12  n 


2i  -  1 1 2 
2  n 


where. 


n  =  sample  size, 

U,  =  F(X(,))  =  CDF(for  normal  distrihu'ion). 


and, 


-^(1)  <  X(2)  <  A'(3)  <  ...  <  Vf„) 
be  n  observations  in  order. 


(2.5) 
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The  C-VM  test  is  also  a  function  of  the  vertical  distance  between  the  hypothesized 
and  the  empirical  distribution  functions  like  the  K-S  test.  But  the  difference  is  that  the  C- 
VM  test  not  only  considers  the  largest  differences  but  also  considers  i!n”  differences  between 
the  two  curves.  Intuitively,  the  C-VM  test  statistic  makes  more  complete  use  of  the  data 
than  K-S  statistic,  but  the  facts  fall  either  to  prove  or  disprove  such  intuition  (5:306).  By 
using  this  fact,  Conover  stated  that  C-VM  test  seemed  intuitively  more  appealing  than 
the  Kolmogorov  test  to  some  people  (5:306). 

Anderson-Darling  brought  the  idea  of  incorporating  a  weight  function  into  K-S  and 
C-VM  statistics  (6:100).  Primarily,  the  Anderson-Darling  statistic  is  based  on  a  weighted 
average  of  the  squared  discrepancy.  Anderson-Darling  test  statistic  is  defined  as  (6:101): 

A-  D  =  -n  -  (l/n)£(2i  -  l)[logt/(i)  +  log(l  -  i/(n_f+1))]  (2.6) 

f=i 


where, 


n  -  sample  size, 

Ui  =  F( A'(ij)  =  CDF(for  normal  distribution), 


-^(1)  <  -X(2)  <  ^(3)  <  •  •  •  <  -Y(n) 

be  n  of-  •  •  vations  in  order.  In  this  expression,  log(a:)  means  loge(:c). 


2.2.2  i'f  .  er  Studies.  The  power  of  a  goodness-of-fit  tea,  is  the  probability  that  the 
test  will  reject  th »  null  hypothesis  (14:473).  The  following  discussion  presents  some  of  the 
modifications  done  for  obtaining  more  powerful  goodness-of-fit  tests. 

Green  and  Hegazy  modified  the  A-D,  C-VM,  and  the  K-S  goodness-of-fit  tests  where 
tin-  parameters  of  the  distribution  were  estimated.  After  the  Monte  Carlo  power  studies, 
they  showed  that  the  modified  A-D  and  C-VM  tests  improved  the  power  of  the  known 
tests  (8:204). 

Porter  modified  the  A-D  and  C-VM  tests  for  the  Pareto  distribution,  lie  used  the  in¬ 
verse  transform  technique  to  generate  the  random  deviates.  Then  he  used  these  deviates  to 


find  the  best  linear  unbiased  estimates  of  the  scale  and  location  parameters.  He  concluded 
that  the  power  of  the  test  improved  when  the  sample  sizes  increased.  He  also  showed  that 
A-D  and  C-VM  tests  were  more  powerful  than  the  Chi-square,  especially  when  the  sample 
data  were  taken  from  the  Weibull,  the  beta,  or  the  normal  distribution  (15:7). 

Ream  modified  K-S,  A-D,  and  C-VM  tests  for  the  normal  distribution  to  test  the 
technique  of  reflecting  data  points  about  the  mean  and  created  tables  of  critical  values  of 
the  modified  K-S,  A-D,  and  C-VM  statistics.  He  used  the  technique  developed  by  Efron 
(1979)  which  was  a  method  for  representing  the  order  statistics  on  a  continuous  spectrum. 
He  applied  this  test  by  plotting  the  values  of  the  order  statistics  and  represented  the  spaces 
between  them  as  piecewise  linear  functions.  He  concluded  that  the  powers  calculated  for 
symmetrical  alternatives  to  the  normal  are  asymptotically  greater  for  the  modified  statistics 
than  for  the  corresponding  unmodified  statistic  values  (17:62). 

Bush  modified  the  A-D  and  C-VM  tests  for  the  Weibull  distribution.  He  presented 
a  Monte  Carlo  method  for  obtaining  the  critical  values  of  the  modified  A-D  and  C-VM 
goodness-of-fit  tests  for  the  three  parameter  Weibull  distribution  when  the  scale  and  loca¬ 
tion  parameters  were  not  specified.  He  also  concluded  that  A-D  and  C-VM  tests  were  not 
very  powerful  when  the  sample  size  was  five.  When  the  sample  size  increased,  he  showed 
that  A-D  and  C-VM  tests  were  more  powerful  than  K-S  test  (3:52). 

Vlviano  modified  the  A-D  and  C-VM  tests  for  the  Gamma  distribution.  He  used 
the  random  deviates  to  estimate  the  maximum  likelihood  scale  and  location  parameters. 
His  power  comparison  showed  that  A-D  goodness-of-fit  test  was  more  powerful  against  the 
lognormal  distribution  (18:47). 

Stephens  concluded  after  his  power  studies  for  normality  and  exponent iality  that 
A-D  statistic  had  better  powers  against  some  alternatives  such  as  the  uniform  and  the 
logistic  distributions  (6:167). 


2.  3  Conclusion 

All  these  research  studies  have  been  done  in  order  to  provide  more  powerful  tests  for 
the  certain  distributions.  They  all  modified  the  test  techniques  or  used  a  new  technique 
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III.  Methodology 


3.  1  Introduction 

This  chapter  initially  discusses  the  thought  processes  associated  with  this  thesis,  and 
then  presents  the  specific  procedures  used  in  generating  the  critical  values  and  the  specific 
procedures  used  in  the  power  tests.  Specifically,  the  discussion  covers  the  topics  of  Monte 
Carlo  simulation  method,  calculation  of  the  critical  values,  the  powers  of  the  goodness-of-fit 
tests,  and  finally  the  details  of  the  computer  programs. 

3.2  Discussion 

3.2.1  Monte  Carlo  Simulation  Method.  The  Monte  Carlo  method  is  used  for  reli¬ 
ability  prediction  when  an  exact  mathematical  model  cannot  be  developed  economically 
or  when  it  becomes  too  complex  to  permit  timely  evaluation.  This  method  involves  the 
determination  of  the  distributions  of  the  various  elements  in  a  system,  selection  of  the 
random  sample  of  each  element,  and  combining  of  these  samples  to  obtain  a  measure  of 
the  system  performance  or  reliability.  The  process  of  random  selection  and  determination 
of  the  system  effects  are  repeated  a  large  number  of  times  and,  each  repetition  results  in 
another  different  estimate  of  the  system  characteristic  that  is  being  measured.  The  law 
of  large  numbers  states  that  ,  as  the  sample  size  increases  ,  the  difference  between  the 
sample  mean  and  the  population  mean  becomes  smaller,  and  the  sample  mean  becomes  an 
increasingly  better  estimate  of  /z.  Then  the  distribution  of  the  sample  becomes  a  better 
representation  of  ihe  distribution  of  the  population  (1:176). 

An  important  characteristic  of  the  Monte  Carlo  method  is  its  reliance  on  computers 
to  simulate  random  processes.  The  Monte  Carlo  approach  is  used  in  this  thesis  to  generate 
the  critical  value  tables  for  the  A-D  and  C-VM  goodness-of-fit  techniques.  The  approach 
is  to  observe  random  variates  chosen  so  that  they  directly  simulate  the  random  processes 
of  the  original  problem.  Ream  tested  the  consistency  of  the  critical  values  with  different 
seeds  each  time  using  150.  300,  500.  1000,  and  5000  samples.  The  only  number  of  samples 
that  yielded  consistent  results  through  thoce  computations  at  all  levels  of  a  was  5000  (17). 
This  thesis  uses  5000  samples  in  the  omputations  as  most  of  the  studies  have  used  in 
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this  area.  Another  reason  why  only  5000  repetitions  were  used  was  to  avoid  using  a  large 
amount  of  computer  time. 

3.2.2  Modified  Test  Statistics.  The  computer  was  used  to  generate  5000  samples 
of  size  n.  Each  group  of  n  normal  deviates  represents  a  simulated  sample  from  the  nor¬ 
mal  distribution.  Then  5000  different  values  of  test  statistics  were  obtained  by  using  the 
modified  A-D  and  C-VM  statistics.  Two  different  modifications  were  used  in  this  thesis: 


1.  Using  the  median  rank  plotting  position 

2.  Using  the  mean  rank  plotting  position 


3.2.2. 1  Using  the  Median  Rank.  There  is  a  lot  of  discussion  in  the  literature 
about  the  choice  of  the  quantile  probabilities.  A  frequently  used  formula  is  given  by  (6: 164): 


Pi  = 


i  -  c 

n  —  2c  +  1 


(3.1) 


where, 

i  =  rank  of  the  order  statistics 
n  =  total  number  of  order  statistics 


c  =  some  constant  satisfying  o  <  c  <  1 


This  research  will  use  c  =  0.3175  since  the  resulting  probabilities  closely  approximate 
medians  of  uniform  (0, 1)  order  statistics.  Using  the  c  value  chosen  ,  the  median  rank  value 
is  given  by: 


i  —  0.3175 
n  +  0.365 


(3.2) 


The  median  rank  equation  given  above  will  be  used  for  the  further  computations. 
Using  the  median  rank  formula  given  above,  the  modified  A-D  and  C-VM  statistics  are 
given  as: 


The  Modified  C  -  V M 


12  n 


+  £ 


t=i 


Ui- 


i-  0.3175 
n  -f  0.365 


(3.3) 


The  Modified  A-D^-n-  +  q  +  ]o^1  ~  ^-+1))]  (3-4) 

The  median  rank  plotting  position  seems  to  be  more  accurate  than  the  other  plotting 
positions  (11:300). 

The  unmodified  C-VM  test  statistic  (2.5)  is  a  function  of  the  values  of  U{  which  is  the 
CDF  value  for  the  normal  distribution  function  U{  is  calculated  at  each  X,  value  in  the 
random  sample,  and  from  this  is  subtracted  the  quantity  (2 i-  l)/2 n,  which  is  the  average 
of  the  empirical  distribution  functions  just  before  and  just  after  the  jump  that  occurs  at 
X ,  ,  that  is,  the  average  of  ( i  -  1  )/n  and  i/n.  That  difference  is  squared  so  that  positive 
differences  do  not  cancel  the  negative  differences,  and  the  results  are  added  together.  After 
this  brief  explanation  of  C-VM  test  statistic,  the  main  purpose  of  the  research  can  be  stated 
again  as,  how  much  more  power  will  be  provided  with  the  new  statistic  using  the  median 
rank  as  the  plotting  position.  The  idea  for  the  modified  A-D  goodness-of-fit  test  is  again 
the  same,  which  is  inserting  the  median  rank  to  the  original  A-D  statistic  < 5:30b). 

3. 2. 2. 2  Using  the  Mean  Rank.  The  other  plotting  position  used  in  this  thesis 
is  the  mean  rank.  The  c  value  in  the  quantile  probability  formula(3.1)  is  equal  to  0  for  the 
mean  rank.  The  mean  rank  formula  is  given  as: 

i 

Mean  Rank  = - -  (3.5) 

n+ 1  ' 

Using  the  mean  rank  formula  given  above,  the  modified  A-D  and  C-VM  statistics 
are  given  as: 

The  Modified  A  -  D  =  -n  -  ^(^^)[logt'r(i)  +  log(  1  -  ,_,+]) )]  (3.0) 

1  =  1 

1  71  7  I2 

The  Modified  C  -  VM  =  —  +  Y  V, - —  (3.7) 

12  tt  “  n  +  1 

i=i 

3. 2. 2. 3  Plotting  Positions.  Given  a  series  of  observations  ordered  from  small¬ 
est  to  largest  (or  vice  versa),  each  event  may  be  assigned  a  plotting  position  which  A  its 
cumulative  probability.  Recurrence  intervals,  which  are  ifciprocah  of  the  rumulnlivr  mob- 
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abilities,  may  ako  be  used  to  define  plotting  positions  (9:i6l5),  The  plotting  positions  are 
used  to  determine  the  percentiles  of  the  distribution  underlying  a  set  of  n  ordered  sample 
values  (9.1615).  The  plotting  position  techniov  involves  using  a  nm  uev  of  discr£  .e  valuta 
of  ih.*  ordered  test  statistics  arid  locating  them  on  a  continuous  spectrum  by  representing 
the  spaces  between  them  as  piecewise  linear  functions  (9:i614).  By  linearly  imerpolatmg 
the  desired  percentiles  between  dk'  mte  values  of  the  test  statistics,  more  accurate  critical 
values  will  be  obtained  (15:1.11). 

'he  cumulative  distribution  function  of  a  sample  size  of  n  is  usually  defined  as  a 
step  function  which  jumps  from  {i  -  1  )/n  to  (i/n)  at  the  ith  order  statistic  of  the  sample. 
When  tile  plotting  position  i/n  is  used,  the  largest  value  cannot  be  plotted.  If  the  plotting 
position  (t’-l)/nis  used,  T  en  the  smallest  value  cannot  bt  slotted.  The  probabilities  0  and 
1  are  off  the  scale  for  rrobabilitv  p«,-2r  constructed  for  tue  normal  distribution  or  for  any 
other  distribution  unlim.v_a  in  extent  (-00  to  oo).  Therefore  statisticians  have  proposed 
diiieren'  plotting  positions,  and  they  have  tried  to  find  the  best  plotting  positions.  Harter 
concluded  that  the  optimum  plotting  position  depended  on  the  use  that  was  to  be  made 
of  the  results  and  also  depended  on  the  underlying  distribution  (9:1615), 

After  his  Monte  Carlo  study  of  plotting  positions,  Harter  found  that  tiro  median 
rank  plot'.ing  position  wai  one  of  the  best  plotting  positions  (10:320).  The  median  rank 
plotting  position,  on  which  this  thesk  places  the  greatest  emphasis,  yields  median  unbiased 
estimates  of  a:,-  for  a  specified  F(xt)  and  of  F(xt)  for  a  specified  x,  (9:1625). 

3.2.3  Calculation  of  the  Critical  Values.,  After  calculating  tiie  5000  statistics  ler 
each  goodness-of-fit  test,  in  order  m  determine  the  critical  values,  it  is  first  necessary  to 
determine  the  critical  regions.  This  critical  region  is  the  set  of  all  points  in  l!  e 
space  which  result  in  the  oecision  to  reject  the  null  hypothesis.  Aker  the  critical  region 
is  determined  ,  then  the  critical  values  can  be  found  with  a  desired  "level  of  significance", 
or  n,  which  is  the  maximum  probability  of  error  in  rejecting  a  tine  null  hypothesis.  If 
hv  is  true,  then  the  maximum  probability  of  error  in  rejecting  II 0  is  a;  and  thereicr*’,  the. 
minimum  probability  of  no  error  in  accepting  //„  is  i  -  a.  In  this  ci.se,  i,!m  value  of  I  -  n 
repicsents  a  certain  percentile  of  the  5000  ordered  test  statistic  vait  os. 


The  technique  used  in  this  thetas  to  determ :  ,e  the  critical  values  is  the  bootstrap 
technique  svhich  was  firs.  *:sed  by  Efron  i?).  Many  staCslldans  have  extensively  uttd  this 
technique  to  determine  the  critical  values. 


The  research  plots  the  n  ordered  test  s' atistir.  values  3i,®3,  . . ii:  along  the  hori¬ 
zontal  axis  and  the  n  plotting  posit1. .-n  values  . .  .  ,yt.  which  are  calculated  with  the 
median  rank  formula  or  mean  rank  formula  alony:  the  vertical  *xis  by  using  t!,!"  tecui:i.;uG, 
These  values  are  assigned  to  positions  2  to  (n  +  1)  on  their  respective  values  on  the  axes. 
There  exists  [0,  t]  interval  on  the  vertical  axis  performed  by  entering  the  endpoints  y$  ~  0 
at  the  position  1  and  yn+1  =  1  at  the  posit'  >n  n  +  2.  Letting  the  order  statistics  be 
represented  by  the  horizontal  <a  ,s  and  lett»*.g  b-  vertical  axis  be  scaled  between  zero  and 
one.  the  statistics  will  be  represent'* on  a  continuoi.'  function.  Then,  the  coritoponding 
endpoints  on  the  horizont?1  ,.xis  arc  found  by  linear  interpolation. 


To  plot  th_  collection  of  u  discrete  values  onto  a  fully  continuous  line  b,  tween  0 
and  3  requires  an  extrapolation  of  the  plotting  axes.  The  first  point  on  the  hoi  'outal 
axis.  x0,  is  computed  by  linearly  extrapolating  from  the  second  and  third  points,  subject 
to  a  non-negativity  restriction.  Extrapolation  is  performed  by  using  the  standard  linear 
slope-intercept  formula  y  =  mx  +  6  to  compute  the  endpoints  a,*o  and  xn+\.  To  find  the 
first  endpoint  on  the  horizontal  axis,  the  slope  is  calculated  by: 


m  = 


U/.-H  ~  Vi) 
(•t’t  +  l  ~  S'!*) 


and  the  intercept  is  found  by: 


S) 


b  -  y,  -  mxi 


(3.9) 


Then  the  lower  endpoint  ,t0  is  found  by: 


Cv  =  (yo  -  b)/m  =  (0  -  l>)  =  -6/m 


■3.10) 


ily  using  the  nonnegativity  rule,  C;>  is  set  to  0.  So  t.nat: 


xq  =  max(Q>-b/m) 


(3.11) 


After  both  endpoints  are  calculated,  n  test  statistic  values  and  their  two  extrapo¬ 
lated  endpoints  are  plotted  on  the  x  -  y  axis.  Then,  the  critical  values  corresponding 
to  the  desired  percentiles  are  found  by  linearly  interpolating  between  the  points  ( x,,y ,) 
and  (a;j+i,yi+i)  using  the  formulas  above;  where  Cp  is  the  critical  value  for  the  desired 
percentile  (15). 

3.2.4  Povy.  0/  ‘ue  6  iodness-of-fit  Tests.  The  goodness  of  a  test  is  measured  by  u 
and  ,5,  the  probabilities  *  r  Type  1  and  Type  II  errors  mentioned  in  Chapter  II,  where  a  is 
chosen  in  advance  and  determines  the  location  of  the  rejection  region.  Basically,  the  power 
of  a  te-*  is  the  probability  that  the  tost  will  reject  the  null  hypothesis. 

Suppose  that  A  is  the  test  statistic  and  RR  is  the  rejection  region  for  a  test  of  a  given 
hypothesis  concerning  the  value  of  a  parameter  9.  Then  the  power  of  the  test,  denoted  by 
power  ((?},  Ik  the  probability  that  the  test  rejects  H0  when  the  actual  parameter  value  is 
9.  That  is, 


power(fl)  =  P  (A  in  RR  when  the  parameter  va'uo  is  9) 


The  power  of  a  test  at  9  —  90  is  equal  to  the  probability  of  rejecting  H0  when  Hc,  is 
true,  where  the  null  hypothesis  is 

H0  ;  9  =  90.  (3.12) 

That  is.  powST(fl0)  =  a.  the  probability  of  a  Type  I  error  (14).. 

The  main  purpose  of  the  power  test  is  to  test  the  thesis  problem  that  the  goodness- 
of-fit  tests  with  the  modified  A-D  and  C-VM  statistics  will  result  in  good ness-of- fit  tests 
having  higher  powers  than  the  powers  of  the  ones  with  the  unmodified  A-D  .md  C-VM 
statistic  values. 

To  execute  the  power  test  ,  random  deviates  from  some  other  alternative  distributions 
of  sample  size  n  described  below  will  be  generated.  Then  the  A-I)  and  C-VM  statistics 


will  he  calculated  with  the  modified  A-D  and  C-VM  statistic  values.  Since  the  critical 
values  are  known  for  the  sample  values  of  n  and  at  the  specified  values  of  a  ,  the  A*l) 
and  C-VM  statistics  calculated  will  be  compared  to  these  respective  values.  This  process 
will  be  repeated  5000  times.  The  number  of  the  times  each  statistic  value  exceeds  the 
respective  critical  value  will  be  counted  for  each  sample  size  of  n.  This  total  number  will 
represent  the  number  of  the  rejections  of  the  null  hypothesis.  This  total  rejection  number 
will  be  divided  by  5000  to  obtain  the  hypothesis  rejection  quotient.  Then  the  power  will 
he  stated  as: 


The  power  at  each  a  —  level  = 


Total  number  of  rejections 
5000 


(3.13) 


After  getting  all  the  powers  at  each  ct-level,  final  conclusions  will  be  derived  for  each 
alternative  distribution. 

The  power  study  is  performed  at  n  =  10,  n  =  20,  n  =  30,  n  =  40,  n  =  50,  n  =  60, 
n  =  70,  and  n  =  80.  By  using  the  numbers  shown  above,  it  will  be  possible  to  compare 
the  powers  of  the  tests  at  different  levels. 


3.2.5  Specific  Procedures.  The  research  effort  will  consist  of  four  phases.  These 
phases  will  be  performed  by  using  the  procedures  and  the  techniques  mentioned  earlier. 
These  phases  are: 


1.  The  A-D  and  C-VM  critical  value  tables  will  be  derived  using  the  known  A-D  and 
C-VM  statistic  values. 

2.  The  A-D  and  C-VM  critical  value  tables  will  be  generated  using  the  modified  A-D 
and  C-VM  statistic  values. 

3.  The  powers  of  the  modified  A-D  and  C-VM  statistics  will  be  determined  by  using 
alternative  distributions.  These  alternative  distributions  are  the  uniform,  the  expo¬ 
nential,  the  double  exponential,  the  Wcibull,  and  the  lognormal  distributions. 

I.  The  conclusions  will  be  derived  about  the  powers  of  the  new  modified  statistics. 
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3.2.5. 1  Generating  the  Critical  Value  Tables.  During  this  phase,  the  critical 
value  tables  will  be  generated  using  the  Monte  Carlo  simulation  of  random  deviates.  The 
following  steps  describe  the  procedure  for  generating  the  critical  value  tables. 

1.,  Random  variates  will  be  generated  for  a  given  sample  size  of  n  from  Normal  distribu¬ 
tion  iV(0, 1).  To  generate  an  array  of  ordered  Ar(0, 1)  random  deviates.the  subroutine 
named  GGNO  will  be  used  from  the  International  Mathematical  and  Statistics  Li¬ 
brary  ( IMSL ). 

2.  The  random  variates  generated  will  be  standardized  by  using  the  following  transfor¬ 
mation: 

2.  =  (3.14) 

a 

These  standardized  data  will  be  used  in  the  further  computations  instead  of  the 
generated  data.Then  the  test  statistics  will  be  computed  from  the  Zl  values  instead 
from  the  original  random  sample. 

3.  The  hypothesized  cumulative  distribution  function  U,  will  then  be  calculated  for 
i  =  1,2,...,  n. 

4.,  The  A-D  and  C-VM  statistics  will  be  calculated  using  the  known  and  the  modified 
statistics. 

5.  Each  of  these  four  steps  described  above  will  be  repeated  5000  times  to  generate  5000 
independent  statistics. 

6.  These  5000  statistics  will  be  ordered  from  smallest  to  largest.  Using  the  plotting 
position  techniciue,  the  80^'^85^h.90^,^95^,l  .and  99i/l  percentiles  of  the  distributions 
of  the  A-D  and  C-VM  statistics  will  be  determined  by  linear  interpolations.  Those 
percentiles  represent  the  .‘20,  .15.  10.  .05.  .01  levels  of  significance.  Then  the  critical 
value  tables  will  be  documented.  These  critical  value  tables  will  be  fully  documented 
to  gu.e  a  cleat  motivation  for  the  future  studies. 

The  flowchart  representation  for  generating  the  critical  values  is  shown  on  the  next 

page. 


3-8 


Figure  3.1.  Generation  of  Critical  Values 
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3. 2.5. 2  Power  Comparison.  In  this  stage  of  the  research  effort,  the  powers  of 
the  modified  A-D  and  C-VM  statistics  will  be  calculated.  As  mentioned  earlier,  the  power 
of  a  test  is  the  probability  that  the  test  will  reject  the  null  hypothesis.  The  null  hypothesis 
is  that  a  sample  of  deviates  follows  a  specified  distribution  where  the  alternative  hypothesis 
is  that  the  sample  of  deviates  follows  some  other  distribution. 

H0  :  Sample  deviates  follow  a  specified  distribution. 

Ha  :  They  follow  some  other  distribution. 

The  steps  for  this  phase  of  the  research  are  as  follows: 

1.  Random  deviates  for  the  selected  distributions  will  be  generated. 

2.  Under  the  given  hypothesis  test  given  above,  the  test  statistics  will  be  compared  to 
the  corresponding  critical  values  in  the  table  for  a  given  level  of  significance  a  .  If 
the  test  statistic  is  equal  to  or  greater  than  the  critical  value  in  the  table,  then  the 
null  hypothesis  ( H0 )  is  rejected. 

3.  The  first  two  steps  will  be  repeated  5000  times,  and  these  rejections  will  be  counted 
to  obtain  the  power  of  each  statistic  value. 

The  flowchart  representation  for  the  power  tests  is  shown  on  the  next  page. 
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Figure  3.2.  Power  Comparison 
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3.3  Conclusion 


This  chapter  represented  the  specific  procedures  used  in  this  research  study.  First 
the  discussion  about  the  Monte  Carlo  simulation  method  was  explained  and  then  the 
techniques  to  generate  the  critical  values  were  presented  with  the  modified  A-D  and  C-VM 
statistics. 

The  next  step  will  be  generating  the  critical  values  and  performing  the  power  tests 
using  the  procedures  explained  in  this  chapter. 
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IV.  Power  Studies 


4.1  Introduction 

This  chapter  presents  the  critical  value  tables  and  the  power  comparison  tables.  The 
methodology  described  in  the  previous  chapter  is  used  to  perform  the  power  comparisons. 
The  complete  tables  of  the  critical  values  are  presented..  Specifically,  the  discussion  covers 
the  topics  of  the  critical  value  tables  and  the  power  comparison  tables. 

4-2  Discussion 

4-2.1  Critical  Value  Tables.  The  complete  tables  of  critical  tables  are  presented  for 
the  following  tests  for  the  sample  sizes  n  =  4  through  n  =  80: 

1.  The  Anderson-Darling  test 

2.  The  Anderson-Darling  test  modified  by  using  the  median  rank 

3.  The  Cramer-von  Mises  test 

4.  The  Cramer-von  Mises  test  modified  by  using  the  median  rank 

The  critical  values  for  the  A-D  and  C-VM  tests  are  generated  by  using  the  known 
A-D  and  C-YM  statistics.  The  critical  values  for  the  modified  A-D  and  C’V-M  tests  are  gen¬ 
erated  by  using  the  new  modified  statistics  given  by  Equations  (3.3).  (3.4).  (3.6).  and(3.7). 
These  four  critical  value  tables  are  included  in  the  study  for  clear  comparisons  between 
the  tests. 
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Table  4.1.  Critical  Values  for  the  Anderson-Darling  Test  (A-D 


Level  of  Significance  (o  = 

Type  I  Error) 

0.20 

0.15 

I 

0.10 

|  0.05  ] 

0.01 

0.4273410 


6  0.436344 


0.4451900 


8  0.4602733 


9  0.4656225 


10  0.4768448 


11  0.4752173 


0.4811802 


3  0.4799352 


14 


15  i  0.483560: 


0.4686368  0.5215520 


0.4730742  0.5266748 


0.5462756 


0.5014591  0.5624942 


0.6195503 


0.8582927 


0.5135141 


0.5213203 


0.5301452 


0.5256196 


0.5274926 


0.5303540 


0.5725465 


0.5830092 


0.5879040 


0.5988379 


0.5938590 


0.596463 


0.6825218 


0.6977583 


0.7005739 


0.7120132 


0.7081370 


0.6977839 


0.4801931 


19  0.4942989 


0.5429153 


0.6148014  0.7455606 


0.5963792  0.7095528 


0.73 


0.6044996  0.7169 


0.9336353 


0.9280258 


0.97-31674 


0.9824843 


0.9436393 


0.9995052 


0.9911392 


0.6116218  I  0. 
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Table  4.2.  Continued  (A-D) 


Level  of  Significance  (o=  Type  I  Error) 


n 

0.20 

0.15 

0.J) 

0.05 

0.01 

22 

0.4930449 

0.5406342 

0.6087568 

0.7233239 

1.004090 

23 

0.4965525 

0.5448771 

0.6276847 

0.7518797 

1.024874 

24 

0.4945974 

0.5435820 

0.6202430 

0.7262725 

0.9730281 

25 

0.4968777 

0.5474281 

0.6153478 

0.7253182 

0.9874505 

26 

0.4870243 

0.5349074 

0.6020432 

0.7138987 

0.9905736 

27 

0.5064116 

0.555  1526 

0.6255779 

0.7485665 

1.035611 

28 

0.5035725 

0.5482607 

0.6243295 

0.7442713 

0.9996597 

29 

0.4967098 

0.5470409 

0.6043720 

0.7331629 

1.018337 

30 

0.5002642 

0.5460176 

0.6153994 

0.7396423 

1.019706 

31 

0.4920015 

0.5412827 

0.6096362 

0.7206756 

0.9634513 

32 

0.4977684 

0.5414925 

0.6100102 

0.7216454 

1.037556 

33 

0.6123351 

0.7156105 

0.9811327 

34 

0.4958267 

0.5410481 

0.6100807 

0.7225074 

1.010662 

35 

0.4973602 

0.5470352 

0.6169090 

0.7334061 

0.9903298 

36 

0.5038891 

0.5516129 

0.6242867 

0.7248192 

0.9966785 

m 

0.4977741 

0.5424100 

0.0118068 

0.7373582 

1.014974 

38 

0.4951401 

0.5434284 

0.6065102 

0.7302265 

0.9693126 

39  0.4948826  0.5474167  0.6189229  0.7319393 


1.035174 


Table  4.3.  Continued  (A-D) 


Level  of  Significance  (a  —  Type  I  Error) 


n 

0.20 

0.15 

0.10 

0.05 

0.01 

40 

0.5036736 

0.5492897 

0.6088142 

0.7162741 

0.9774227 

41 

0.4935474 

0.5456238 

0.9670739 

42 

0.4987202 

0.5480480 

0.6173229 

0.7275543 

1.032360 

43 

0.5078812 

0.5625439 

0.6357059 

0.7474669 

1.015957 

44 

0.5057907 

0.5523701 

0.6249388 

0.7439404 

1.030191 

45 

0.5075741 

0.5600471 

0.6272904 

0.7511843 

1.049332 

46 

0.5009937 

0.5491180 

0.6143340 

0.7353191 

1.079781 

0.5439053 

0.6071224 

0.7375965 

1.043659 

48 

0.5051956 

0.5512581 

0.6226959 

0.7505187 

1.048113 

49 

0.4975300 

0.5474415 

0.6134128 

0.7330531 

0.9863949 

50 

0.5064582 

0.5549145 

0.6306038 

0.7509441 

1.017109 

51 

0.5016975 

0.5513611 

0.6221276 

0.7502614 

1.011376 

52 

0.5054054 

0.5616665 

0.6354789 

0.7535743 

1.013035 

53 

0.4991112 

0.5465317 

0.6146659 

0.7347736 

0.9884982 

54 

0.4980010 

— ■ 

0.7294196 

0.9858552 

55 

0.5054150 

jlllllll 

0.6312294 

0.7499961 

1 .044233 

56 

0.5054855 

0.5514259 

0.6198158 

0.7405357 

1.001810 

57 

0.5003051 

0.5516874 

0.6244812 

0.7361297 

1.019163 
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Table  4.4.  Continued  (A-D) 


L 

evel  of  Significance  (a  — 

Type  I  Errol 

r) 

0.20 

0.15 

0.10 

0.05 

0.01 

0.5041008 

0.5561218 

0.6305560 

0.7591591 

1.017233 

0.4955597 

0.5473194 

0.6220034 

0.7387219 

0.9823399 

0.5019493 

0.5448074 

0.6112021 

0.7257367 

1.023222 

0.4970607 

0.5500755 

0.6146068 

0.7381287 

1.029202 

0.5038795 

0.5488949 

0.6174755 

0.7389317 

0.9899883 

0.5008945 

0.5522652 

0.6225338 

0.7408447 

0.9984367 

0.5083809 

0.5651475 

0.6351047 

0.7575722 

1.035576 

0.5034600 

0.5512009 

0.6210823 

0.7369995 

1.019421 

0.5056152 

0.5513268 

0.6228599 

0.7395401 

0.9776191 

0.4956703 

0.5476418 

0.6135138 

0.7234917 

1.002739 

0.4987144 

0.5423241 

0.6143417 

0.7315444 

0.9953613 

0.4891967 

0.5374489 

0.6148720 

0.7375908 

1.016052 

0.5006447 

0.5449486 

0.6086006 

0.7174299 

0.9944189 

0.5036469 

0.5571595 

0.6310691 

0.7516361 

0.9956664 

0.5042152 

0.5534744 

0.6187972 

0.7528877 

1.098568 

0.4993744 

0.5456963 

0.6249349 

0.7475891 

0.9832656 

0.4938010 

0.5401573 

0.6125678 

0.7294769 

1.027989 

0.5070686 

0.5569229 

0.6184120 

0.7290077 

1.064887 
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Table  4.5.  Continued  (A-D) 


■ 

Level  of  Significance  (a-Type  I  Error) 

_ 19 

n 

0.20 

0.15 

0.10 

0,05 

0.01 

76 

0.5072517 

0.5548745 

0.6221275 

0.7382354 

0.9811439 

77 

0.5040436 

0.5481758 

0.6155281 

0.7394637 

1.029260 

m 

0.5019684 

0.5506668 

0.6074256 

0.7241516 

1.010888 

79 

0.4939728 

0.5408135 

0.6081811 

0.7273674 

1.041717 

80 

0.5535432 

0.6213608 

0.7532310 

1.046329 

Table  4.6.  Critical  Values  for  the  Modified  Anderson-Darling  Test  (Mod.  A-D) 


Level  of  Significance  (a=Type  I  Error) 


n  0.20  0.15  0.10  0.05 


IS  I  1.595449 


0.05 

0.01 

1.748723 

1.901822 

1.786559 

2.001122 

1.786873 

2.006685 

1.801773 

2.026566 

1.827742 

2.073187 

1.832235 

2.107482 

1.821748 

2.080344 

1.835227 

2.068199 

1.832728 

2.093671 

1.841662 

2.111019 

1.834993 

2.112786 

1.821638 

. 

2.075947 

1.837376 

2.129101 

1.866905 

2.117235 

1.829279 

2.056677 

1.848537 

2.097478 

1.836084 

2.122495 

1.838660 

2.080507 

Table  4.7.  Continued  (Mod.  A-D) 


Level  of  Significance  (q= 

Type  I  Error) 

0.20 

0.15 

0.10 

0.05 

0.01 

1.605027 

1.654015 

1.723168 

1.837928 

2.123885 

1.607399 

1.656692 

1.740698 

1.866060 

2.144170 

1.654651  I  1.731568  1.840148  2.089002 


1.838591  2.104431 


1.606775  1.657915 


1.596288  1.645237  j  1.714891 


1.614830  1.665790  1.735625  1.862446  2.148955 


1.612313  1.657729  1.735475  1.854534  2.111488 


1.604225  1.655478  1.714090  1.842946  2.132516 


i 

1.607267  1.654267  1.724511  1.850140  |  2.132860 


1.598879  1.648886  1.717329  1.830383  2.074648 


—mum i—w 


1.604487 


1.602903  1.654169  1.719116 


1.601793  1.648222  1.716944  1.830750  j  2.120502 


2.100450 


1.608988  1.657530  1.729572  1.832800  !  2.105381 


— I— —I 


1.602711 


1.599356 


1.712374  1.836834  I  2.077459 


1.599217  1.651529  ;  1.724192  |  1.837242  i  2.143903  II 
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Table  4.8.  Continued  (Mod.  A-D) 


Level  of  Significance  (g=  Type  I  Error) 


0.20  0.15  !  0.10 


0.05  0.01 


1.608160  1.654381  1.713497  1.822373  2.084799 


1.597586  1.649262  1.714718  j  1.832047  2.075401 


1.602928  1.652935  1.722067  1.833937  2.139005 


1.611481  1.666656  1.740896  1.854076  2.122475 


1.608974  1.657139  1.729172  1.848476  2.138180 


1.610071  1.663834  1.731770  1.854791  2.154166 


1.603903  1.652592  1.718451  1.838369  2.184988 


1.604656  1.654627  1.728684  1.855408  2.123127 


1.602488  1.654079  1.726776  1.846361  2.114164 


1.597328  1.644581  1.706415  1.833174  2.111530 


1.606836  1.661249  1.734949  1.861311  2.148796 


.GU9284  1.659309  1.738390  1.839491  2.131019 


1.603733  1.654411  1.722675  1.845199  2.132654 


1.603451  1.652371  1.713387  1.824596  2.170734 


1.598221  1.649824  1.720283  1.833208  2.110764 


1.604191  1.658361  1.733399  1.844929  2.126796 


1.600681  1.653149  1.724302  1.841839  2.134912 


1.604868 


1.835598  2.124434 


Table  4.9.  Continued  (Mod.  A-Dj 


Level  of  Significance  (a  —  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

58 

1.608841 

1.658838 

1.720072 

1.842489 

2.122024 

59 

1.601927 

1.653996 

2.123772 

1 

1.604246 

1.850096 

2.099027 

61 

1.601652 

1.649261 

1.721050 

1.842003 

2.114895 

62 

1.609081 

1.664414 

1.726728 

1.855917 

2.137005 

63 

1.606510 

1.658436 

EH 

1.851661 

64 

1.595528 

1.643341 

1.708275 

1.821423 

2.082424 

65 

1.614925 

1.665764 

1.738296 

1.848522 

2.120224 

66 

1.603245 

1.655537 

1.723713 

1.853168 

2.155743 

|  67 

1.605503 

1.652283 

1.727871 

1.854118 

2.106010 

68 

1.601856 

1.646454 

EH 

1.821701 

2.072914 

69 

1.602928 

1.650131 

1.721455 

1.826721 

|  70 

1.601944 

1.646496 

1.711170 

1.825813 

I  *  1 

EH 

1.644791 

1.713329 

1.848534 

II 

2.133702  | 

72 

1.604134 

1.651684 

1.718487 

73 

1.593395 

1.643261 

1.709183 

1.825314 

- r 

2.140033 

74 

1.602291 

1 .652550 

1.724125 

1.857033 

2.137829 

i 

75 

1.594906 

1.645405 

1.721660 

1.841076 

2.113717 

Table  4.10.  Continued  (Mod.  A-D) 


Level  of  Significance  (a= 

Type  I  Error) 

11 

» 

0.15 

0.10 

0.05 

0.01 

1.608738 

1.661160 

1.740814 

1.863881 

2.192662 

— 

1  TT 

i _ 

1.599689 

m 

1.725956 

1.831131 

2.113838 

!-« 

i 

1.598564 

1.648224 

1.718864 

1.836502 

2.126775 

i 

79 

■9 

1.719360 

80 

1.602745 

1.650009 

2.109204 
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Table  4.11.  Critical  Values  for 


■ 

Level  of  Signifi 

n 

0.20 

0.15 

4 

6.06S717889 

0.075482644 

5 

0.074497014 

0.81921719 

6 

0.073810749 

0.081494279 

B 

0.075224698 

0.082806356 

8 

0.077101640 

0.085106254 

9 

0.078383930 

10 

0.079117730 

0.087914445 

11 

0.078922927 

0.087871805 

12 

0.09329304 

0.088333301 

13 

0.078199118 

0.088464923 

14 

0.0/898i8/4 

0.087575160 

15 

0.079056472 

0.087549478 

1 16 

0.080250062 

0.088360421 

B 

0.079889290 

0.090032421 

18 

0.086529292 

19 

0.081109054 

0.089803994 

20 

21 

0.079958476 

0.089077838 

the  Cramer-vt  n  Mises  Test  (C-VM) 


cance  (q=  Type  I  Error) 


0.10 

0.05 

0.01 

0.085986301 

0.1027864 

0.1329755 

0.092780255 

0.1088244 

0.1519448 

■ 

0.1578436 

0.097887017 

0.1198000 

0.1676036 

0.1204440 

0.1719634 

0.1196247 

0.1692951 

0.1000359 

0.1211092 

0.1659990 

0.1021110 

0.1221171 

0.17x6552 

0.1737344 

0.099740401 

0.1191476 

0.1690551 

1— 

0.1211501 

0.1750606 

i 

|  0.1039800 

0.1289154 

0.1754416 

B— 1 

0.1663492 

B 

0.1242585 

0.1726213 

0.1816628 

0.1013554 

DEB 

Table  4.12.  Continued  (C-VM) 


n  0.20 


Level  of  Significance  (o=  Type  I  Error) 


0.088679813  0.1018423  |  0.1238366  0.1773855 


23  0.080085248  0.090274885 


24  0.07939031  0.089777172 


Em 


hhhiiimih«uih«bihii 


0.1747164 


25  0.080698572  0.090523042  0.1030444  I  0.1230224  0.1699734 


26  0.078952312  0.087606333 


0.1706200 


0.081380449  0.090854980  0.1030964  0.1283479  0.1768174 


28  I  0.081017911  0.089316107 


0.1728324 


29  0.080269307  0.089592382  0.1017790  0.1247838  0.1799154 


30  0.080438256  0.089790776  0.1020064  0.1237940  0.1782377 


31  i  0.079447605  0.088063426 


32  |  0.079861253  0.089245029 


ItlilTt 


0.1808827 


0.089789651  0.1020840  I  0.1211182  I  0.1760548 


0.088616826 


IflHl 


0.089403555  0.1023602  ;  0.1261002  j  0.1755391 


0.089757077 


Line 


37  I  0.081238300  0.088818684  0.1010038  j  0.1264580  1  0.1787453 


38  i  0.080053426  0.088507704  |  0.1006704  i  0.1224907  I  0.1633722 


39  |  0.080003247  0.088^44657  0.1019481 
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Table  4.13.  Continued  (C-VM) 


■ 

Level  of  Significance  (a=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

40 

0.080477841 

0.089769222 

0.1011164 

0.1196890 

0.1695593 

41 

0.080171429 

0.088017531 

0.1005671 

0.1244502 

0.1731949 

42 

0.080914006 

0.089308478 

0.1018946 

0.1228404 

0.1791724 

43 

0.081907399 

0.091332987 

0.1750055 

44 

0.081682011 

0.091009520 

0.1021398 

0.1265333 

0.1841971 

45 

0.081103258 

0.090990670 

0.1033571 

0.1269290 

0.1839232 

46 

0.080394760 

0.089128010 

0.1025386 

0.1233494 

0.1827025 

47 

0.079413734 

0.088695556 

0.1002844 

0.1236093 

0.1751411 

48 

0.080372520 

0.090312369 

0.1031824 

0.1251452 

0.1786935 

49 

0.080017462 

0.089598760 

0.1016602 

0.1224438 

0.1748529 

50 

Q.08001880 1 

0.089651220 

0.1030667 

0.1253582 

0.1759703 

51 

0.080380484 

0.089901276 

0.1032750 

0.1250169 

0.1754528 

52 

0.080479883 

0.090384185 

0.1038498 

0.1822425 

53 

0.088755876 

| 

0.1016509  !  0.1230755 

0.1731220 

55 

0.090576157 

0.1043266 

0.1260005 

0.1849225  j 

56 

0.090332493 

QEB 

57 

0.080314606 

0.089454904 

0,1027690 

0.1240410 

0.1762333  i 
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Table  4.14,  Continued  (C-VM) 


Table  4.15.  Continued  (C-YM) 


Level  of  Significance  {o=  Type  I  Error? 


76  0.081169255  0.090258129  0.1024938  0.1219029  0.1646494 


0.081446439  0.090075724  0.1035862  0.1237362  0.1 


0.081172861  0.090091988  |  0.1031678 


0.080688916  0.089598104  0. 


7  !  0.17744 


80 

0.079321444 

.089231446 

1 

0.1024222 

0.1257416 

0.1837831 

-I- lb 


Table  4.16.  Critical  Values  for  the  Modified  Cramer-von  Mises  Test  (Mod.  C-VM) 


Level  of  Significance  (a~  Type  I  Error) 


n  0.20 


0.1171122  0.1258052  0.1399408  0.1631716 


0.1239169  0.1363130  0.1598520 


6  0.1068020  0.1162525  0.1293700  0.1531270 


0.1032833  0.1128287  0.1268466  0.1502074 


8  0.1016296  0.1121755  0.1274291  0.1523345 


9  0.1002520  0.1098706  0.1242078  0.1462136 


10  0.096721902  0.1067272  0.1205297 


11  0.094881043  U.1041535  0.1188708  0.1445517 


12  0.094152793  0.1041171  0.1179808  0.1418343 


13  0.094214745  0.1040925  0.1186687 


U  0.094684064  0.1043528  0.1173616  0.1403915 


15  0.094134524  0.1043178  0.1183724  0.1416568 


16  0.092405163  0.1024797  0  1 176013  0.1412600 


17  0.090981118  0.1004940  0.1148321  0.1395562 


0.01 


0.2042842 


0.2155906 


0.2002948 


0.1996575 


0.2109509 


0.2049321 


0.2029956 


0.2027317 


8 

0.089875177 

0.098631233 

9 

0.090443790 

0.099732995 

0.1973863 


0.1953897 


0.1948619 


0.2043033 


0.2003421 


20  |  0.089838564  0.099070512  0.1129879  j  0.1362366 


21  0.089684360  0.099597573 


0.1969647 


0.1878682 


Table  4.17.  Continued  (Mod.  C-VM) 


Level  of  Significance  (a=  Type  I  Error) 


n 

0.20 

0.15 

0.10 

0.05 

0.01 

22 

■■ 

0.099156566 

0.1116720 

0.1353266 

0.1956280 

23 

0.0892  o3 

0.097931020 

0.1120457 

0.1355183 

0.1845801 

24 

0.090250760 

0.098927520 

0.1116665 

0.1345997 

0.1894391 

25 

0.086477816 

0.096268803 

0.1097936 

0.1330310 

0.1883914 

26 

0.089361295 

0.099474996 

0.1116120 

0.1350340 

0.1839367 

27 

0.087897420 

0.096669361 

0.1099983 

0.1325055 

0.1888126 

28 

0.087610364 

0.097102761 

0.1115020 

0.1350892 

0.1910306 

29 

0.086007208 

0.096212253 

0.1096980 

0.1325047 

0.1878328 

30 

0.086151421 

0.095770031 

0.1086233 

0.1290127 

0.1805344 

31 

0.085716859 

0.096077800 

0.1094933 

0.1336658 

0.1866961 

32 

0.087037697 

0.097366460 

0.1106475 

0.1361363 

0.1887586 

33 

0.086482324 

0.1297646 

0.1838285 

34 

0.086857319 

0.096330233 

0.1084818 

0.1334283 

0.1921022 

35 

0.086611502 

0.096350551 

0.1082830 

0.1293105 

0.1893613 

36 

0.1291684 

0.1835041 

37 

0.085253775 

0.094642401 

0.1082638 

0.1305203 

0.1826830 

38 

0.084169395 

0.093680479 

0.1060074 

0.1283366 

0.1774915 

39 

0.084604226 

0.1302742 

0.18.30927 
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Table  4.18.  Continued  (Mod.  C-VM) 


Level  of  Significance  (o—  Type  I  Error) 


n 

0.20 

0.15 

0.10 

0.05 

0.01 

40 

0.087177113 

0.096347749 

0.1103945 

0.1337443 

0.1834177 

41 

0.085038245 

0.094490454 

0.1059541 

0.1316031 

0.1858505 

42 

0.084178187 

0.093363374 

0.1072612 

0.1291613 

0.1772463 

43 

0.086647049 

0.096142627 

0.1099741 

0.1314089 

0.1859473 

44 

0.086165667 

0.096147276 

0.1089791 

0.1321802 

0.1807189 

45 

0.086863168 

0.097464323 

0.1116782 

0.1368531 

0.1855542 

46 

0.084294438 

0.094070919 

0.1069317 

0.1279515 

0.1944696 

47 

0.083825737 

0.093185492 

0.1074729 

0.1312952 

0.1849413 

48 

0.084716879 

0.093679935 

0.1077310 

0.1310238 

0.1805953 

49 

0.084885985 

0.093539461 

0.1057767 

0.1298060 

0.1897912 

50 

0.084612489 

0.093931831 

0.1081517 

0.1312491 

0.1883270 

51 

0.085586265 

0.095265456 

0.1085853 

0.1319058 

0.1854176 

52 

0.085252240 

0.094202630 

0.1075557 

0.1286628 

0.1816324 

53 

0.082118019 

0.092360727 

0.1055925 

54 

0.084748983 

0.093850747 

0.1066248 

0.1306041 

0.1898644 

55 

0.085000128 

0.093799673 

0.1080838 

0.1314944 

0.1822819 

56 

0.084181368 

0.093914039 

0.1338640 

57 

0.082907148 

0.091027714 

i 

0.1043118 

0.1271668 

0.1732585 

Table  4.19.  Continued  (Mod.  C-YM) 


Level  of  Significance  (o=  Type  I  Error) 

0.20 

0.15 

0.10 

0.05 

0.01 

0.082885757 

0.092641443 

0.1053403 

0.1266784 

0.1769804 

0.083405033 

0.093811668 

0.1067430 

0.1295628 

0.1896901 

0.083464295 

0.093149699 

0.1059194 

0.1299687 

0.1804228 

0.081543289 

0.090840355 

0.1036204 

0.1264776 

0.1685492 

0.084958017 

0.095521927 

0.1087770 

0.1315241 

0.1841160 

0.082068309 

0.092810810 

0.1069565 

0.1284048 

0.1864124 

0.082610682 

0.091001593 

0.1046564 

0.1266130 

0.1744784 

0.082847416 

0.092799S72 

0.1055391 

0.1267786 

0.1766667 

0.081755854 

0.090612829 

0.1042599 

0.1255016 

0.1724503 

0.082927428 

0.093032897 

0.1063054 

0.1294967 

0.1764327 

0.083560646 

0.093030259 

0.1060542 

0.1276196 

0.1837067 

0.083318591 

0.093147844 

0.1069267 

0.1281051 

0.1858886 

0.085421734 

0.094470188 

0.1080142 

0.1290280 

0.1830784 

0.083139651 

0.093213074 

0.1076766 

0.1280835 

0.1768681 

0.083336644 

0.092905864 

0.1057336 

0,1265798 

0.1828810 

0.084840626 

0.093278311 

0.1057599 

0.1302512 

0.1799054 

0.085315138 

0.093913630 

0.1075264 

0.1314406 

0.1812564 

0.084850922 

0.094617531 

0.1083276 

0.1322037 

0.1838400 
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Table  4.20.  Continued  (Mod.  C-VM) 


■ 

Level  of  Significance  (a=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

76 

0.082804367 

0.091524869 

0.1058075 

0.1247443 

0.1755374 

1 7 

0.082122885 

0.091925353 

0.1047585 

0.1266247 

0.1827777 

78 

0.084189937 

0.093079522 

0.1064301 

0.1291159 

0.1803580 

79 

0.084658489 

0.094539940 

0.1080402 

0.1278412 

0.1756680 

80 

0.083445959 

0.093562253 

0.1065563 

0.1313581 

0.1870539 
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4-2. '2  Power  Comparison  Tables.  This  section  contains  two  sets  of  power  compari¬ 
son  tables.  The  first  set  of  tables  compares  the  powers  of  the  tests  modified  by  using  the 
median  rank  to  the  powers  of  the  known  A-D  and  CY-\1  goodness-of-fit  tests.  The  second 
set  of  tables  compares  the  powers  of  the  tests  modified  by  using  the  mean  rank  to  the  known 
A-D  and  CY-M  goodness-of-fit  tests.  These  two  sets  of  power  comparisons  are  presented 
for  each  alternative  distribution  to  display  the  powers  of  the  modified  goodness-of-fit  test1 
clearly. 

In  addition  to  these  power  tables,  the  visual  representations  for  the  power  improve¬ 
ments  are  given  for  each  distribution  with  each  plotting  position  used.  These  figures  show 
the  power  improvements  as  the  sample  size  doubles  for  the  significance  level  of  o  =  0.05 
which  is  the  most  commonly  used  (4). 
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Table  4.21.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Uniform  Distri¬ 
bution  Using  the  Median  Rank 


Level  of  Significance  (a—  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01  j 

10  T 

0.2806000 

0.2090000 

0.1510000 

10  { 

0.2828000 

0.2128000 

0.1516000 

0.0836000 

_ 

0.0146000 

20  f 

0.4804000 

0.3962000 

0.2904000 

20  i 

0.4846000 

30  f 

0.6394000 

0.5648000 

0.4526000 

0.2906000 

0.0828000 

30  { 

0.6436000 

0.5676000 

0.4564000 

40  T 

0.7516000 

0.6914000 

0.6036000 

40  t 

0.7542000 

0.6930000 

0.6082000 

50  f 

0.8568000 

0.5614000 

i 

0.2682000 

50  1 

0.8002000 

0.7148000 

0.5562000 

0.2480000 

60  f 

0.8278000 

| 

0.7072000  |  0.3774000 

60  i 

0.9154000 

i 

70  f  |  0.9548000 

0.9354000 

0.9002000 

i 

0.8140000  I  0.5272000 

BUB 

80  t  j  0.9794000 

0.8502000 

;  0.5952000 

warn 

0.9384000 

0.8740000 

0.6316000 

f Unmodified  A-D  Test 
{Modified  A-D  Test 


Table  4.22.  Power  Comparison  for  the  Cramer- von  Mises  Test  with  the  Uniform  Distri¬ 
bution  Using  the  Median  Rank 


■ 

Level  of  Significance  (o=  Type  I  Error) 

n 

0.20 

n 

0.15 

0.10 

0.05 

0.01 

10  f 

0.2676000 

0.2018000 

0.1422000 

0.0741999 

0.1180000 

EB 

0.2570000 

0.1202000 

0.0557999 

0.0820000 

:  22  f 

0.4262000 

0.3462000 

0.2496000 

0.1452000 

0.0233999 

HJ 

0.3630000 

0.2980000 

0.2100000 

0.1156000 

0.0179999 

m 

0.5560000 

0.4688000 

0.3714000 

0.2342000 

0.0575999 

5 

0.5180000 

0.4344000 

0.3376000 

0.2164000 

0.0575999 

!  40  t 

0.6552000 

0.5780000 

0.4888000 

0.3578000 

0.1268000 

mm 

Efli 

■Mil 

0.6006000 

0.5272000 

0.4228000 

0.2796000 

0.0916000 

!  50  f 

| 

0.7594000 

0.6924000 

0.5896000 

U.4266000 

0.1730000 

50  t 

0.7288000 

0.6568000 

0.5422000 

0.3814000 

0.1288000 

60  f 

i 

0.8388000 

0.7820000 

0.7056000 

0.5468000 

0.2878000 

Q1 

0.8080000 

0.7470000 

0.6604000 

0.4984000 

0.2272000 

j  70  f 

0.8876000 

0.8454000 

0.7660000 

0.6344000 

0.3174000 

m 

0.8638000 

0.8154000 

0.7296000 

0.5984000 

0.2874000 

m 

0.9324000 

0.8316000 

0.7082000 

0.3736000 

m 

0.9100000 

0.8674000 

i 

0.8044000  !  0.6668000 

0.3434000 

fUiunodified  C'-VM  Test 
•  Modified  C'-VM  Test 
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Table  4.23.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Weibull  Distribu- 


Table  4.24.  Power  Comparison  for  the  Cramer- von  Mises  Test  with  the  Weibull  Distri 


bution  Using  the  Median  Rank 


Level  of  Significance  (o=  Type  I  Error) 


0.20 

0.15 

0.10 

0.4380000 

0.3768000 

0.3062000 

20  t  0.6922000  [  0.6322000  0.5518000  |  0.4458000  |  0.2362000 


l 

0.7568000  0.7092000  i  0.6364000  0.5192000  0.3016000 


30  t  [  0.8086000  i  0.7700000 


0.3988000 


30  |  j  0.8674000 


40  t 


•  jHjl 


00  0.7800000  j  0.7002000  0.5064000 


i 

0.8272000  !  0.7596000  0.5714000 


fiillllllglliPWOilllltllRiWiPIlUl 


RIiPlIiItfcMiliriJ 


40  i  I  0.9288000  0.9034000  0.8646000  0.7862000  0.6176000 


50 


!  |  i 

50  t  I  0.9622000  j  0.9482000  0.9244000  I  0.8676000  0.7124000 

!  1  | 


0.9728000  |  0.9624000  j  0.9424000  ■  0.9040000  !  0.7926000 


0.9820000  0.9744000  |  0.9602000  |  0.9248000  1  0.8154000 


0.9846000 


I 

I 

I 


0.9538000 

0.8756000 

0.9614000 

0.8816000 

0.9952000  j  0.9924000  I  0.9872000  j  0.9732000  .  0.9116000 


flnmodified  C-YM  Test 
iModified  C-YM  Test 
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Figure  4.3.  Anderson-Darling  Power  Test  for  the  Weibull  Distribution  with  the  Median 
Rank 


Figure  4.4..  Cramer-von  Mises  Test  for  the  Weibull  Distribution  with  the  Median  Rank 


Table  4.25.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Exponential  Di 
tribution  Using  the  Median  Rank 


Level  of  Significance  (q=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  t 

0.6662000 

0.6046000 

0.5332000 

0.4236000 

0.2228000 

m 

0.6660000 

0.6058000 

0.5338000 

0.4240000 

0.2240000 

20  f 

0.9280000 

0.9026000 

0.8590000 

0.7866000 

0.5778000 

20  i 

0.9290000 

0.9024000 

0.8588000 

0.7538000 

0.5780000 

30  f 

0.9868000 

0.9822000 

0.9682000 

0.9366000 

0.8128000 

•l-l- 

o 

CO 

mnuu 

0.9822000 

0.9684000 

0.9366000 

0.8120000 

40  f 

0.9984000 

0.9892000 

0.9332000 

40  | 

0.9892000 

mm 

m 

0.9996000 

0.9964000 

0.9822000 

m 

0.9996000 

0.9996000 

0.9792000 

0.9996000 

ISHIlB 

^80 

0.9992000 

0.9956000 

n 

0.9996000 

0.9996000 

0.9996000 

0.9992000 

0.9964000 

70  t 

1.000000 

1.000000 

m 

IBB 

1.000000 

1.000000 

1.000000 

0.9990000 

m 

1.000000 

1.000000 

1.000000 

BBS 

0.9994000 

iB!™|!m!bwimwim 

f Unmodified  A-D  Test 
iModified  A-D  Test 


Table  4.26.  Power  Comparison  for  the  Cramer-vori  Mises  Test  with  the  Exponential  Dis 
tributiori  Using  the  Median  Rank 


Level  of  Significance  (o=  Type  I  Error 


0.20 

0.15 

0.10 

0.6402000 

0.5746000 

0  5044000 

0.7354000 

0.6808000 

0.6026000 

0.8982000 

0.8666000 

0.8200000 

IMHftCiXfitl 


0-93 


I 

0.9762000  I  0. 

1 


I 

0.9874000  0.9800000  |  0.9666000  0-9392000 


TjgjgiimgjgW 


0-795*000 


0.9466000  0.9010000 


0.9958000 


0.9976000 


0.9990000 


0.9880000  0-975400.  j  0.9112000 


IQOliliyTiTiTilgtlt 


0.9814000  0.931 


0.99<3800Q  0-9910000 


iiM 


0.9994000  j  0.9992000  0.9982000  0-9948000  0.9716000 


0.9992000  0.9988000  0.9976000  0. 


0.9996000  0.9994000  i  0.9992000  0.9986000  0.9930000 


.000000  I  1.000000  j  1.000000 


.000000  I  1. 000000  j  1 .000000 


1.000000  I  1.000000  0.9996000  !  0.99x200 

1  S 


1.000000  !  l.OOOOun  i  l.OOOnOO  |  0.999X000  |  0.99XX000 


Table  4.27.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Lognormal  Dis 
tribution  Using  the  Median  Rank 


Level  of  Significance  (a=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  t 

0.4562000 

0.3890000 

0.3210000 

0.2356000 

0.1046000 

10  + 

0.4550000 

0.3908000 

0.3220000 

0.2354000 

0.1048000 

20  f 

0.7246000 

0.6708000 

0.5950000 

0.4808000 

0.2860000 

20  t 

0.7252000 

0.6694000 

0.5948000 

0.4420000 

0.2850000 

30  t 

0.8392000 

0.8050000 

0.7536000 

0.6518000 

0.4524000 

30  t 

0.8402000 

0.8050000 

0.7530000 

0.6512000 

0.4514000 

40  t 

0.9266000 

0.9016000 

0.8700000 

0.8056000 

0.1656000 

40  + 

0.9268000 

0.9016000 

0.8706000 

0.8044000 

0.6348000 

50  i 

0.9590000 

0.9466000 

0.9238000 

0.8712000 

0.7346000 

50  i 

0.9604000 

0.9458000 

0.9236000 

0.8682000 

0.7234000 

60  t 

0.9826000 

0.975800 

0.9462000 

0.9352000 

0.8244000 

60  | 

0.9818000 

0.9742000 

0.9612000 

0.9304000 

0.8402000 

70  f 

0.9912000 

0.9878000 

0.9810000 

0.9660000 

0.9034000 

70  { 

0.9910000 

0.9878000 

0.9804000 

0.9650000 

0.8998000 

80  t 

0.9976000 

0.9950000 

0.9912000 

0.9802000 

0.9290000 

80  i 

0.9976000 

0.9952000 

0.9920000 

0.9822000 

0.9380000 

fUnmodified  A-D  Test 
IModified  A-D  Test 


Table  4.28.  Power  Comparison  for  the  Cramer-von  Kises  Test  with  the  Lognormal  Dis¬ 
tribution  Using  the  Median  Rank 


■ 

Level  of  Significance  (a=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  f 

0.4450000 

0.3772000 

0.3080000 

0.2112000 

0.0726000 

0.5960000 

0.5294000 

0.4468000 

0.3300000 

0.1488000 

20  f 

0.6922000 

0.6322000 

0.5518000 

0.4458000 

0.2362000 

m 

0.7568000 

0.7092000 

0.6364000 

0.5192000 

0.3016000 

30  f 

0.8220000 

0.7780000 

0.7196000 

0.6242000 

0.3896000 

30  t 

0.8516000 

0.8186000 

0.7630000 

0.6710000 

0.4468000 

40  f 

0.9060000 

0.87*  ; 

0.8278000 

0.7504000 

0.5286000 

++ 

o 

— 

0.9194000 

0.8970000 

0.8562000 

0.7780000 

0.5710000 

50  f 

0.9490000 

0.9296000 

0.8978000 

0.8422000 

0.6568000 

m 

0.9530000 

0.9392000 

0.9154000 

0.8564000 

0.6822000 

E9 

0.9738000 

0.9630000 

0.9420000 

0.9064000 

0.7590000 

m 

0.9770000 

0.9682000 

0.9514000 

0.9146000 

0.7794000 

70  f 

0.9862000 

0.9800000 

0.9676000 

0.9438000 

0.8314000 

m 

0.9866000 

0.9838000 

0.9716000 

0.9468000 

0.8434000 

m 

0.9930000 

0.9894000 

0.9836000 

0.9664000 

0.8850000 

eg 

0.9936000 

0.9906000 

0.9846000 

0.9684000 

0.8948000 

fUnmodified  C-YM  Test 
I  Modified  C-YM  Test 
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10  20  30  40  50  60  70  80 

Sample  Size 


Figure  4.7.  Anderson-Darling  Power  Test  for  the  Lognormal  Distribution  with  the  Me¬ 
dian  Rank 


10  i  30  40  50  60  70  80 

Sample  Size 


Figure  4.8.  Cramer-von  Mises  Power  Test  for  the  Lognormal  Distribution  with  the  Me¬ 
dian  Rank 


Table  4.29.  Power  Comparison  for  the  Anderson- Darling  Test  with  the  Double  Exponen 
tial  Distribution  Using  the  Median  Rank 


■ 

Level  of  Significance  (a=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  f 

0.T352000 

0.6T98000 

0.6130000 

0.5002000 

0.2878000 

lot 

0.7358000 

0.6810000 

0.6132000 

0.5004000 

0.2888000 

20  f 

0.9588000 

0.9440000 

0.9144000 

0.8618000 

0.6896000 

20  t 

0.9592000 

0.9440000 

0.9150000 

0.8354000 

0.6894000 

30  f 

0.9950000 

0.9918000 

0.9856000 

0.9666000 

0.8968000 

30  i 

0.9950000 

0.9918000 

0.9858000 

0.9664000 

0.8960000 

40  t 

0.9996000 

0.9996000 

0.9992000 

0.9964000 

0.7180000 

40  1 

0.9996000 

0.9996000 

0.9992000 

0.9964000 

0.9774000 

50  t 

1.000000 

1.000000 

1.000000 

0.9992000 

0.9950000 

50  t 

1.000000 

1.000000 

1.000000 

0.9992000 

0.9938000 

60  t 

1.000000 

1.000000 

1.000000 

1.000000 

0.9992000 

60  t 

1.000000 

1.000000 

1.000000 

1.000000 

0.9996000 

TO  t 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

TO  t 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

80  f 

| 

1.006000 

1.000000 

1.000000 

1.000000 

1.000000 

80  1 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

fUnmodified  A-D  Test. 
JModified  A-D  Test 
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Table  4.30.  Power  Comparison  for  the  Cranier-von  Mises  Test  with  the  Double  Exponen 
tial  Distribution  Using  the  Median  Rank 


Level  of  Significance  (a=  Type  1  Error) 


0.9890000 


tiX&MHIHil  ■HJtfiilMiHJ 


1.000000 


i. 000000 


1.000000 


0.05 

0.01 

0.4662000 

0.2624000 

0.5678000 

0.3x38000 

0.8258000 

0.6158000 

0.8670000 

0.6880000 

0.9464000 

0.8556000 

0.9690000 

0.9020000 

0.9894000 

0.9598000 

0.9912000 

0.9688000 

0.9966000 

0.9854000 

0.9984000 

0.9888000 

0.9996000 

0.9978000 

0.9998000 

0.9986000 

1.000000 

0.9992000 

1.000000 

0.9996000 

1.000000 

0.9998000 

1.000000 

1.000000 

t Unmodified  C'-VM  Test 
i.Modified  C-VM  Test 
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Figure  4.9.  Anderson-Darling  Power  Test  for  the  Double  Exponential  Distribution  with 
the  Median  Rank 


Figure  4.10.  Cramer- von  Mises  Power  Test  for  the  Double  Exponential  Distribution  with 
the  Median  Rank 
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Table  4.31.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Uniform  Distri 
bution  Using  the  Mean  Rank 


■ 

Level  cf  Significance  (o  = 

i 

Type  I  Error)  j 

n 

0.20 

0.05 

0.01 

m 

0.2806000 

0.2090000 

0.1510000 

0.0821999 

0.0146000 

10  i 

0.0152000 

20  f 

0.4804000 

0.3962000 

0.0374000 

20  1 

0.5103000 

0.4012000 

0.3124000 

0.1891000 

0.0472100 

m 

0.6394000 

0.5648000 

0.4526000 

0.2906000 

0.0828000 

30  1 

0.6488000 

0.5770000 

0.4660000 

0.3024000 

0.0869999 

40  f 

0.7516000 

0.6914000 

0.6036000 

0.4612000 

0.0799999 

40  1 

0.7610000 

0.6962000 

0.1878000 

50  f 

0.8568000 

0.8018000 

0.7126000 

0.5614000 

0.2682000 

50  1 

0.3008000 

60  f 

i 

0.8278000 

0.7072000 

0.3774000 

m 

0.9100000 

0.6766000 

0.3626000 

70  f 

0.5272000 

70  1 

0.9552000 

0.7948000 

0.4902000 

1 

j  80  f 

0.8502000 

0.5952000 

■ 

0.9796000 

0.8626000 

0.5714000 

fUmnodified  A-D  Test 
iModified  A-D  Test 
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Table  4.32.  Power  Comparison  for  the  Cramer- von  Mises  Test  with  the  Uniform  Distri 
bution  Using  the  Mean  Rank 


■ 

Level  of  Significance  (o  = 

Type  I  Error) 

■ 

0.20 

0.15 

OUO 

0.05 

0.01 

10  f 

0.2676000 

0.2018000 

0.1422000 

0.0741000 

0.0118000 

10  I 

0.3312000 

0.2688000 

0.1904000 

0.1056000 

0.0223999 

20  f 

0.4262000 

0.3462000 

0.2496000 

0.1452000 

0.0233999 

20  t 

0.5166000 

0.4310000 

0.3304000 

0.2078000 

0.0524000 

30f 

0.5560000 

0.4688000 

0.3714000 

0.2342000 

0.0575999 

30  | 

0.6498000 

0.5706000 

0.4640000 

0.3172000 

0.0913999 

40  f 

0.6552000 

0.5780000 

0.4888000 

0.3578000 

0.1268000 

0.6692000 


0.6924000 


0.7588000 


0.7820000 


0.8394000 


0.8454000 


0.8880000 


0.8944000 


0.9206000 


0.5752000 


0.5896000 


0.6726000 


0.7056000 


0.7708000 


0.7660000 


0.8224000 


0.8316000 


0.8710000 


0.4266000  0.1730000 


0.5138000  I  0.2420000 

; 


0.5468000  0.2878000 


0.6260000  |  0.3556000 


0.6984000  I  0.3814000 


0.7082000  !  0.3736000 
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Table  4.33.,  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Weibull  Distribu¬ 
tion  Using  the  Mean  Rank 


|i 

Level  of  Significance  (a-  Type  I  Error)  . 

n 

0.20 

0.15 

0.10 

0.01  ! 

1 

m 

0.2020000 

0.1518000 

0.1082000 

0.0572000 

0.0114000 

! 

m 

0.2020000 

0.1538000 

0.1106000 

m 

0.2544000 

0.1966000 

0.1382000 

0.0700000 

0.0162000 

1  20  i 

0.2610000 

0.1988000 

0.0107200 

m 

0.2706000 

0.2110000 

0.1486000 

0.0806000 

0.0210000 

30  | 

0.2748000 

0.2132000 

0.1520000 

0.0816000 

0.0208004  | 

! 

Efl 

0.2976000 

0.2358000 

0.1774000 

0.1074200 

0.0297000  I 

42» 

° 

0.2996000 

0.2352000 

0.0254000 

i 

i  50  f 

ii  _ 

0.3266000 

0.2646000 

0.1860000 

0.1068000 

0.0284000 

IS 

0.3452000 

0.2876000 

!  60  t 

0.3712000 

0.3072000 

0.2336000 

0.1396000 

0.0342000  | 

IS 

0.3522000 

0.2846000 

0.2064000 

Bill 

0.4006000 

i\ 

0.0511999  | 

!  *0  i 

0.3964000 

0.3310000 

|  80  f 

0.4296000 

0.3548000 

0.2778000 

0.1600000 

0.0478000  I 

:  so  j 

• 

0.4246000 

0.3536000 

0.2722000  s  0.1686000 

I 

0.0397999  !' 

| 
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f  Unmodified  A-D  Test 
^Modified  A-D  Test 


Table  4.34.-  Power  Comparison  for  the  Cramer- von  Mises  Test  with  the  Weibull  Distri¬ 


bution  Using  the  Mean  Rank 


Level  of  Significance  (o  = 

Type  I  Error) 

0.20 

0.15 

0.10 

0.05 

0.01 

0.1986000 

0.1558000 

0.1092000 

0.0562000 

0.0112000 

0.0124000 


0.2482000  0.1922000  0.1290000  0.0692000  0.0124000 


frlUUilKT 


0.2630000 


0.2740000  0.2204000 


0.0140000 


0-0766000  0.0189999 


0.0206000 


0.2828000  0.2154000  0.1638000  0.1022000  0.0250000 


0.0280000 


!  ! 

0.3140000  0.2472000  I  0-1752000  |  0.1004000  0.0255999 


j  50  l 

0.3378000 

0.1076000 

0.0286000  j 

I  60  f 

0.3474000 

0.1180000 

0.0392000  i 

0-3728000 


0.2286000  j  0-12-38000  0.0388000 


0.1302000  0.0399999 


0.0392000 


0.3978000  i  0.3208000  j  0.2158000  !  0.1428000  0.0357999 


0.3976000  !  0.3234000  j  0.2482000  i  0.1516000  0.0399999 
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10  20  30  40  50  60  70  80 

Sample  Size 


Figure  4.13.  Anderson-Darling  Power  Test  for  the  Weibull  Distribution  with  the  Mean 
Rank 


Fieure  4.14.  Cramer-von  Mises  Power  Test  for  the  VVeibuli  Distribution  with  the  Mean 


Table  4.35.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Exponential  Dis 
tribution  Using  the  Mean  Rank 


Level  of  Significance  fo=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  f 

0.6662000 

0.6046000 

0.5332000 

0.4236000 

0.2228000  I| 

10  | 

0.6656000 

0.6062000 

0.5358000 

0.4260000 

0.2234000  || 

20  f 

0.9280000 

0.9026000 

0.8590000 

0.7866000 

Q.o « 1 8000  1 

1 

|  20  | 

0.930000 

0.9028000 

0.8610000 

0.7930000 

0.5800000  I 

30  f 

0.9868000 

0.9822000 

0.9682000 

0.9366000 

0.8128000  j 

! 

30? 

0.9868000 

0.9822000 

0.9694000 

0.9368000 

ft 

0.8120000  | 

40  f 

0.9984000 

0.9974000 

0.9958000 

0.9892000 

0.9332000  j 

40  i 

0.9984000 

0.9974000 

0.9948000 

0.9842000 

0.9500000  || 

50  f 

0.9996000 

0.9996000 

0.9992000 

0.9964000 

| 

0.9822000  j 

50  i 

0.9996000 

0.9996000 

0.9994000 

0.9376000 

0.9854000  I 

60  I 

0.9996000 

0.9996000 

0.9996000 

0.99^2000 

0.9956000  | 

j  60  i 

0.9996000 

0.9996000 

0.9996000 

0.9990000 

0.9950000  I 

I  TO  f 

1.000000 

1.000000 

1.000000 

1.000000 

0.9990000  | 

o 

1.000000 

1.000000 

1.000000 

1.000000 

0.9990000  I 

|  80  1  .  1 .000006 

1.000000 

1.000000 

1.000000 

0.9994000  I 

|  80  ? 

'*0000 

1.000000 

1.000000 

1 .000000 

0.9994000  | 

f Unmodified  A-D  Test 
JModined  A-D  Test 
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Table  4.36.  Power  Comparison  for  the  Cramer-von  Mises  Test  with  the  Exponential  Dis¬ 
tribution  Using  the  Mean  Rank 


■ 

Level  of  Significance  (o=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  f 

0.6402000 

0.5746000 

0.5044000 

0.3898000 

0.2012000 

10  + 

0.6352000 

0.5734000 

0.4938000 

0.3766000 

0.1952000 

20  f 

0.8982000 

0.8200000 

0.7424000 

0.4960000 

20  | 

0.8340000 

0.8624000 

0.8134000 

0.73840CC 

0.5072000 

30  | 

0.9762000 

0.9466000 

0.9010000 

0.7490000 

30  t 

0.9760000 

0.9652000 

0.9446000 

0.9004000 

0.7406000 

40  f 

0.9880000 

0.9754000 

0.9112000 

40  | 

0.9960000 

0.9932000 

0.9872000 

0.9740000 

0.9084000 

50  j 

0.9990000 

0.9968000 

0.9910000 

0.9640000 

50  t 

0.997000C 

0.9906000 

0.9640000 

0.9992000 

0.9988000 

0.9976000 

0.9916000 

60  t 

0.9992000 

0.9992000 

0.9988000 

0.9978000 

0.9908000 

1.000000 

1.000000 

1.000000 

1.000000 

0.9958000 

70! 

1.000000 

1.000000 

1.000000 

1.000000 

0.9958000 

80  j 

1.000000 

1.000000 

1.000000 

0.9996000 

0.9982000 

80  ! 

1.000000 

1.000000 

i.ooodfeo 

0,9996000 

0.9984000 

fUnmodified  C-VM  Test 
IModified  C-VM  Test 


10  20  30  40  50  60  70  80 

Sample  Size 


Figure  4.15.  Anderson- Darling  Power  Test  for  the  Exponential  Distribution  with  the 
Mean  Rank 


Power 


Figure  4.16.  Cramer- von  Mises  Test  for  the  Exponential  Distribution  with  the  Mean 


Table  4.37.  Power  Comparison  for  the  Anderson- Darling  Test  with  the  Double  Exponen 
tial  Distribution  Using  the  Mean  Rank 


Level  of  Significance  (a=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

iot 

0.7352000 

0.6798000 

0.6130000 

0.5002000 

0.2878000 

10  { 

0.7352000 

0.6814000 

0.6138000 

0.5012000 

0.2884000 

20  f 

0.9588000 

0.9440000 

0.9144000 

0.8618000 

0.6896000 

20  { 

0.9588000 

0.9442000 

0.9146000 

0.8700000 

0.6898000 

30  j 

0.9950000 

0.9918000 

0.9856000 

0.9666000 

0.8968000 

30  t 

0.9950000 

0.9918000 

0.9860000 

0.9664000 

0.8956000 

40  f 

0.9996000 

0.9996000 

0.9992000 

0.9964000 

0.7180000 

40  % 

0.9996000 

0.9996000 

0.9992000 

0.9946000 

0.9748000 

50  t 

1.000000 

1.000000 

1.000000 

0.9992000 

0.9950000 

50  i 

1.000000 

1.000000 

1.000000 

0.9996000 

0.9962000 

60  t 

1.000000 

1.000000 

1.000000 

1.000000 

0.9992000 

60  t 

1.00000U 

1 ,000000 

1.000000 

1.000000 

0.9996000 

70  t 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

O 

++ 

1.000000 

1.000000 

1.000000 

1.000000 

0.9998000 

80  j 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

80  l 

1.000000 

1.000000 

1.000000 

1.000000 
_ _ 

1.000000 

{Unmodified  A-D  Test 
{Modified  A-D  Test 


Table  4.38.  Power  Comparison  for  the  Cramer- von  Mises  Test  with  the  Double  Exponen¬ 


tial  Distribution  Using  the  Mean  Rank 


Level  of  Significance  (a=  Type  I  Error) 


0.9890000 


HI&BE  OHIUli  lOrfcliUMi] 


0.05 

0.01 

0.4662000 

0.2624000 

0.4502000 

0.2582000 

0.8258000 

0.6158000 

0.8174000 

0.6158000 

0.9464000 

0.8556000 

0.9444000 

0.8486000 

0.9894000 

0.9598000 

0.9886000 

0.9582000 

0.9966000 

0.9854000 

0.9966000 

0.9860000 

0.9996000 

0.9978000 

0.9996000 

0.9968000 

1.000000 

0.9992000 

1.000000 

0.9992000 

1.000000 

0.9998000 

1.000000 

1.000000 

tUnmodified  C-VM  Test 
iModified  C-VM  Test 
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Figure  4.17.  Anderson-Darling  Power  Test  for  the  Double  Exponential  Distribution  with 
the  Mean  Rank 


Figure  4.18.  Cramer-von  Mises  Power  Test  for  the  Double  Exponential  Distribution  with 
the  Mean  Rank 
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Table  4.39.  Power  Comparison  for  the  Anderson-Darling  Test  with  the  Lognormal  Dis 
tribution  Using  the  Mean  Ra\k 


■ 

Level  of  Significance  (o=  Type  I  Error) 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  t 

0.4562000 

0.3890000 

0.1046000 

10| 

0.4534000 

0.3896000 

0.1048000 

20  f 

0.7246000 

0.6708000 

0.5950000 

0.4808000 

0.2860000 

m 

0.7300000 

0.6711000 

0.5990000 

0.4812000 

0.2861000 

j  30  f 

0.8392000 

0.8050000 

0.7536000 

0.6518000 

0.4524000 

0.8040000 

0.7534000 

0.6518000 

0.4482000 

40  T 

0.9266000 

0.9016000 

0.8700000 

0.8056000 

0.6156000 

40  | 

0.9268000 

0.9006000 

0.8636000 

0.7760000 

0.6192000 

50  f 

0.9590000 

0.9466000 

0.9238000 

0.8712000 

0.7346000 

m 

0.9516000 

0.9318000 

0.8864000 

0.7490000 

0.9758000 

0.8244000 

60  1 

0.9726000 

0.9582000 

0.9254000 

0.8128000 

0.9878000 

0.9810000 

0.9660000 

0.9034000 

i 

nn 

0.9906000 

0.9870000 

0.9792000 

0.9602000 

0.8910000  j 

m 

0.9912000 

0.9802000 

0.9290000 

■ 

0.9974000 

0.9946000 

0.9222000 

fUnmodified  A-D  Test 
tModified  Test 
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Table  4.40.  Power  Comparison  for  the  Cramer-von  Mises  Test  with  the  Lognormai  Dis 
tribution  Using  the  Mean  Rank 


■ 

Level  of  Significance  (a=  Type  I  Error)  J 

n 

0.20 

0.15 

0.10 

0.05 

0.01 

10  f 

0.4380000 

0.3768000 

0.3062000 

0.2192000 

0.0939999 

m 

0.4282000 

0.3642000 

0.2912000 

0.2004000 

0.0851999 

20  f 

0.6922000 

0.6322000 

0.5518000 

0.4458000 

0.2362000 

m 

0.6790000 

0.6124000 

0.5408000 

0.4288000 

0.2294000 

30  f 

0.8086000 

0.7700000 

0.7156000 

0.6146000 

0.3988000 

o 

CO 

0.8010000 

0.7586000 

0.7036000 

0.6050000 

0.3814000 

40  f 

0.8996000 

0.8694000 

0.8272000 

0.7596000 

0.5714000 

40  i 

0.8936000 

0.8666000 

0.8206000 

0.7430000 

0.5574000 

50  f 

0.9456000 

0.9272000 

0.8926000 

0.8312000 

0.6748000 

m 

0.9438000 

0.9202000 

0.8872000 

0.8220000 

0.6678000 

n 

0.9728000 

0.9624000 

0.9424000 

0.9040000 

0.7926000 

!  60  i 

i 

0.9718000 

0.9596000 

0.940*000 

0.8992000 

0. 1 1 78000 

IQ 

0.9846000 

0.9786000 

0.9650000 

0.9364000 

0.8398000 

m 

0.9758000 

0.9632000 

0.9334000 

0.8294000 

m 

0.9924000 

0.9888000 

0.9820000 

0.9614000 

0.8816000 

i  80  | 

0.9914000 

0.9884000 

0.9798000 

0.9602000 

0.8824000 

tUnmodified  C-VM  Test 
1  Modified  C-VM  Test 
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Figure  4.19.  Anderson-Darling  Power  Test  for  the  Lognormal  Distribution  with  the  Mean 
Rank 


1 

0.9 
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0.7 

Power  0.6 
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Figure  4.20.  Cramer- von  Mises  Power  Test  for  the  Lognormal  Distribution  with  the  Me¬ 
dian  Rank 
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Sample  Size 


4-3  Conclusion 


This  chapter  is  mainly  a  collection  of  tables.  The  critical  values  are  presented  for  the 
known  and  the  modified  tests  for  the  sample  sizes  n  =  4  through  n  =  80  at  the  significance 
levels  of  g  =  .20.  .15.  .10.  .05.  .01.  The  power  comparison  tables  contain  two  sets  of  tables. 
These  compare  the  powers  of  the  modified  tests  to  the  known  A-D  and  CV*M  tests.  The 
results  of  the  power  comparisons  will  be  given  in  the  next  chapter. 


V.  Results 


5.1  Introduction 

This  chapter  discusses  the  results  of  the  power  comparisons  presented  in  the  previous 
chapter  and  the  computer  programs  used.  Specifically,  the  discussion  covers  the  topics  of 
the  computer  programs,  the  test  of  the  random  number  generator,  the  validation  and 
verification,  and  the  analysis  of  the  power  comparisons. 

5.2  Discussion 

5.2.1  Computer  Programs.  There  are  two  main  programs  used  in  the  study.  One 
of  them  is  for  calculating  the  critical  values,  and  the  other  is  for  performing  the  power 
comparisons.  Both  programs  are  written  in  FORTRAN'  77  and  are  run  in  AFITNET 
Mainframe.  In  addition  to  these  programs,  some  subroutines  are  taken  from  the  Interna¬ 
tional  Mathematical  and  Statistics  Library(IMSL).  Both  of  these  programs  are  presented 
in  the  Appendices.  The  detailed  comments  are  given  to  help  the  readers  understand  the 
logic  of  the  programs  clearly.  The  flow  diagrams  for  the  computer  programs  are  given  in 
Chapter  3. 

5.2.2  Test  of  the  Random  Number  Generator.  To  generate  the  critical  value  ta¬ 
bles  precisely,  it  is  essential  to  use  a  reliable  random  number  generator.  Forty  numbers 
were  generated  by  using  the  subroutine  GGNO  from  the  International  Mathematical  and 
Statistics  Library  (IMSL).  GGNO  generates  an  array  of  ordered  N(0.1)  random  deviates. 
The  method  of  normal  probability  plotting  was  used  to  investigate  the  goodness-of-fit  of 
these  deviates  generated  to  ‘he  normal  distribution.  After  the  deviates  were  plotted,  a  line 
was  drawn  ’‘by  eye"  through  these  points  (6:28).  Two  points  on  the  plot  corresponding 
to  around  the  10£"  percentile  and  another  around  the  90tft  percentile  were  located  and 
connected  to  each  other.  This  line  showed  a  very  good  fit. 

5.2.8  Validation  and  Verification.  The  validation  of  the  computer  programs  is  per¬ 
formed  by  comparing  the  critical  values  to  the  other  published  critical  values.  The  critical 
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values  for  the  known  A-D  and  C-VM  tests  were  compared  to  the  ones  obtained  by  Ream 
(17).  The  values  were  the  same  for  the  unmodified  goodness-of-fit  tests. 

The  purpose  of  the  verification  is  to  ensure  that  the  concepts  and  the  equations 
derived  in  the  study  are  accurately  reflected  by  the  computer  programs..  To  achieve  this, 
every  line  of  the  computer  programs  was  checked  many  times  by  the  author. 

0.2.4  The  Analysis  of  Power  Comparisons.  Using  the  median  rank  plotting  po¬ 
sition  improved  the  power  of  the  C-VM  goodness-of-fit  test  approximately  10 %  for  the 
lognormal,  the  exponential,  the  Weibull.  and  the  double  exponential. 

Using  the  mean  rank  plotting  position  also  improved  the  power  of  the  C-VM  goodness- 
of-fit  test  only  for  the  uniform  and  the  Weibull  distributions:  while  it  reduced  the  powers 
for  the  other  three  distributions.  The  improvement  for  the  uniform  distribution  was  9%: 
while  the  improvement  was  2%  for  the  Weibull  distribution. 

Using  the  median  rank  plotting  position  improved  the  power  of  the  A-D  goodness- 
of-fit  test  for  each  of  the  distributions  used.  The  improvements  were  approximately  at  the 
level  of  0.01%. 

Using  the  mean  rank  plotting  position  also  improved  the  power  of  the  A-D  goodness- 
of-^t  test  for  the  uniform,  the  double  exponential,  and  the  Weibull  distributions.  The 
power  remained  almost  the  same  for  the  exponential  and  the  lognormal  distributions.  The 
improvements  for  the  uniform,  the  double  exponential,  and  the  Weibull  distributions  were 
0.01%. 

The  following  table  shows  the  approximate  power  improvements  of  the  modified  tests 
over  the  unmodified  tests. 
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Table  5.1.  Power  Improvements  of  Modified  Tests  over  Unmodified  Tests 


A-D 

C*  YM 

i  Median 

Mean 

i 

Median  I  Mean 

Uniform 

.01% 

.01% 

5-0% 

9.0% 

Wei  bull 

.01% 

.01% 

10.0% 

2.0% 

Exponential 

.01% 

No  change 

10.0% 

Reduced 

Lognormal 

.01% 

No  change 

o 

© 

Reduced 

Double  Exp. 

.01% 

.01% 

10.0% 

Reduced 

5.3  Conclusion 

Both  plotting  positions  improved  the  powers  of  the  known  A-D  and  C*VM  goodness- 
of-fit  tests.  The  median  rank  plotting  position  seemed  to  have  better  powers  than  the 
mean  rank  plotting  position  did. 


17.  Conclusions  and  Recommendations 


6.1  Conclusions 

Based  on  the  results  obtained  in  this  thesis,  the  following  conclusions  are  noted: 

1.  The  power  comparison  study  based  on  the  alternative  distributions  showed  that  both 
plotting  positions  increased  the  power  of  the  Anderson-Darling  goodness-of-fit  test  for 
almost  aU  the  distributions.  The  power  remained  the  same  only  for  the  exponential 
arid  the  lognormal  distributions  when  the  mean  rank  plotting  position  was  used. 
The  improvement  for  the  C-VM  goodness*of-fit  test  was  considerably  higher  for  the 
uniform  and  the  WeibuU  distributions  when  either  the  median  or  the  mean  rank 
plotting  position  was  used  and.  for  the  lognormal,  the  exponential,  the  WeibuU.  and 
the  double  exponential  distributions  when  the  median  rank  plotting  position  was 
used.  The  median  rank  plotting  position  showed  a  higher  power  than  the  mean  rank 
plotting  position. 

2.  The  critical  value  tables  for  both  tests  were  presented.  The  tables  can  be  used  to 
test  whether  a  random  sample  of  data  foUows  the  normal  distribution. 

6.2  Recommendations 

Based  on  the  observations  made  throughout  the  thesis  study,  the  following  research 
areas  are  proposed  for  further  studies: 

1.  Apply  the  goodness-of-fit  techniques  developed  in  this  thesis  for  other  distribution 
functions. 

2.  Further  modify  the  goodness-of-fit  techniques  developed  in  this  thesis  to  allow  for 
unknown  values  of  the  parameters  of  the  hypothesized  normal  distribution. 
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Appendix  A.  Compnt'r  Programs 


C ***«****♦***************»*♦**»#***♦******»-<*= ********************* 


C* 

C* 

C* 

C* 

C* 

C* 


THIS  PROGRAM  DETERMINES  THE  CRITICAL  /ALUES  WITH  * 

THE  MODIFIED  ANDERSON-DARLING  AND  CRAHER-VOM  MISES  STATISTICS  * 
VALUES  FOR  THE  NORMAL  DISTRIBUTION.  * 


* 


C*****************************«*x********************************** 

c 

c 

c  **  SOME  COMMENTS  ABOUT  THE  VARIABLES  USED  ** 

C 

C  S  *  Number  of  repetitions 

C  M  =  Sample  size 

C  AIDER(*)  =  A-D  Statistic  values 
C  CRIT(*)  *  Critical  values  at  (*) 

C  SEED  =  Input/Qutput  double  precision  triable  assigned  an  integer 
C  value  in  the  exclusive  range 

C  i,N  (GGNO)  *  “I"  is  the  first  order  statistic  to  be  generated, 

C  "N"  is  the  last  order  statistic  to  be  generated 

C  XX(*)  =  Output  vector  of  length  [H+l-IFIRSTj  containing  the  order 
C  statistics  in  an  ascending  order 

C  IER  =  Error  parameter  output 

C  ZZ  =  Input  value  at  which  function  is  to  evaluated 

C  YY  =  Output  probability  that  a  random  variable  having  a  normal 

C  (0,1)  distribution  will  be  less  than  or  eguai  to  ZZ. 

INTEGER  S,J,M,K,SS 

REAL  XX(IOO) ,ZZ,TT,ANDER(0:5010) ,DARL 
REAL  CRITSO , CRIT6S , CRIT90 , CRIT95 , CR1TS9 
C 

C  Enter  the  double  precision  seed 
C 

DOUBLE  PRECISION  SEED 
SEED=  469857936. DO 
C 

C  Number  of  repetition  is  5000 
C 

S=50Q0 


C  enter  the  number  sf  sample  size 
r 


N=8G 


DO  11  J=1,S 


C 

C  Call  GGNO  to  generate  set  of  order  statistics 
C  from  the  normal  distribution 

C 

CALL  GGN0(SEED,1,N,N,XX,ERR) 

C 

C  Call  STANDAR  to  standardize  the  random  deviates 
C  by  Equation(3. 12) .  The  following  computations  will  use 
C  these  stadardize  values. 

C 

CALL  STANDAR(XX.N) 

C 

C  Call  VSRTA  to  sort  the  arrays  by  algebraic  values 
C 

CALL  VSRTA (XX, N) 

DO  22  K=1  ,N 
ZZ-XX(K) 

C 

C  Call  MDNOR  to  determine  the  CDF  values  of  the  normal 
C  distribution  for  the  standardized  data 

C 

CALL  MDNOR(ZZ.TT) 

XX(K)*TT 

22  CONTINUE 

C 

C  Call  ANDERSON  to  find  the  A-D  statistic  values 
C 

C  **Call  CRAMER  to  find  the  CV-M  statistic  values 
C 

CALL  ANDERSON (N, XX, DARL) 

C 

C  »*  CALL  CRAMER (N, XX, DARL )  ** 

C 

ANDER(J)=DARL 
11  CONTINUE 

ANDER(0)=0.0 
SS=  S+l 
C 

C  Call  VSRTA  to  sort  the  statistics  to  determine  the  percentiles 
C  of  the  array  in  order  to  find  the  critical  values 
C 

CALL  VSRTA (ANDER.SS) 

C 

C  Call  EXTRA  to  extrapolate  the  data  for  the  points  at  "0", 
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and  "n+l".  Output  vector  is  an  array  of  length  "n+2". 


CALL  EXTRA (S, ANDER) 

Call  VALUES  to  determine  the  critical  values 
CALL  VALUES (ANDER , CR1T80 , CRIT85 , CRIT90 , CRIT95 , CRIT99 , S) 
Print  the  critical  values 

WRITE ( * , * )  CRIT80, CRIT85 , CRIT90 , CRIT95 , CRIT99 
END 


Using  the  technique  explained  in  Chapter  3,  find  the  critical  values 

SUBROUTINE  CAN(Y1,Y2,D1,D2,Y,RES) 

REAL  M»B,Y1,D1,Y2,02,Y,RES 
M=  (Y2-Y1)/(D2-D1) 

B  =  Yl  -  (M*D1) 

RES  *  (Y-B)/M 
END 


The  following  subroutine  calculates  the  A-D  statistic  values 


SUBROUTINE  ANDERSON (K,R, ANDER) 

INTEGER  K,I 

REAL  R(*), TOTAL, XX, YY, ANDER, ZZ,RR 
TOTAL  *  0.0 
DO  11  1*1, K 
RR  *  R(I) 

IF  (RR  .  LT.  0.0  .AND.  RR.EQ.O.O)  THEN 
RR  =  .0001 
END  IF 

ZZ  *  1.0  -  R(K-I+1) 

IF (ZZ . LE . 0 . 0) THEN 
ZZ  *  .0001 

ENDIF 

XX  *  LOG(RR) 

YY  *  LOG(ZZ) 

TOTAL  =  ( ( (2 . 0*REAL ( I ) ) -0 . 635) * (XX+YY) ) +TOTAL 
CONTINUE 
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ANDER  =  -  REAL(K)  -  (Cl .0/ (REAL 00-0. 365)  )*TGTAL) 
END 


C 

C  The  following  subroutine  calculates  the  CV-M  statistic  values 
C 

SUBROUTINE  CRAfc£R(N » R , ANDER) 

INTEGER  I,N 

REAL  R(*) .TOTAL, ANDER, MEDIAN 
TOTAL  =0.0 
DO  22  1*1, N 

MEDIAN  =  (REAL(I)-O. 3175)/ (REAL(H)-O. 365) 

TOTAL  *  TOTAL  +  (CR(I)  -  MEDIAN)  *  (R(I)  -  MEDIAN)) 

22  CONTINUE 

ANDER  *  (1.0  /  (12.0  *  REAL(N)))  +  TOTAL 
END 

C 

C  Standardize  all  the  data 

C 

SUBROUTINE  STANDAR(X,N) 

INTEGER  I ,N 

REAL  X(*),XSUM.XBAR.S,XQUT 
XSUM  =  0.0 
XOUT  *  0.0 
DO  100  I  *  1,N 

XSUM  =  XSUM  +  X(I) 

100  CONTINUE 
C 

C  Compute  the  mean 

C 

XBAR=XSUM/N 
DO  200  I  =1 ,N 

XOUT  =  XOUT+((X(I)-XBAR)*(X(I)-XBAR)) 

200  CONTINUE 

S  =  SqRT(XOUT/(N-l)) 

DO  300  I  *  l.N 
C 

C  Compute  the  standardized  data 

C 

X(I)  =  (X(I)-XBAR)/S 
300  CONTINUE 
END 
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SUBROUTINE  EXTRA(N.D) 

INTEGER  N,NO,Ni 

REAL  Yi ,Y2,D(0:*) ,D1 ,D2,ZZ 

Yl  -  0.5/N 

Y2  *  i.o/N 

D1  =  D(i) 

D2  =  D(2) 

CALL  CAN(Y1 ,Y2,D1 ,D2,0.0,ZZ) 
IF(ZZ.GE.O.O)  THEN 
D(O)  =  ZZ 
ELSE 

D(O)  *  0.0 
EMDIF 

Yl  =  CREAL(N)  -1.5)/M 
Y2  =  (REAL CM)  -0.5)/N 
NO  *  N-l 
D1  *  D(NO) 

D2  *  DCN) 

CALL  CAN(Yl,Y2,Di,D2»1.0,ZZ) 
N1  *  N+l 
D(N1)  *  ZZ 
END 


C 

C  The  following  subroutine  determines  the  percentiles 

C  and  finds  the  critical  values  by  evoking 
C  the  subroutine  CAN 

C 


SUBROUTINE  VALUES  (D ,  CRIT80 ,  CRIT85 ,  CRIT90 ,  CRUSH ,  GRITS* ,  M's 
INTEGER  I.N.NN 

REAL  D(0:*) ,Y(0:6000) ,C80,C90,C95,C99,C8S, 

1  Y79 ,D79 , Y81 , D81 ,DIF90 ,  Y89, Y91 , D89 ,D91 , DIF95 . DIF80 , 

1  Y94,Y96 , D94 , D96 ,DIF99 , Y98 , Y100 ,D98 ,D100 ,DIF85 , 

1  Y84 , D84 , Y86 , D86 , CRIT8S , CRIT80 , CRIT90 , CRIT9S , CRIT99 

DO  100  1*1, N 

Y(I)  *  (REAL(I)  -  0.5)/REAL(N) 

100  CONTINUE 
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Y(0)  =  0,0 
HN  -  N  +  1 

ycnn)  -  i.o 

C80  ■  1000, 0 
085  -  1000,0 
090  =  1000,0 
095  *  1000.0 
099  =  1000.0 
DO  200  I  =  NN,0,«1 

IF  (Y(I) .LE.0.75)  GO  TO  300 
IF  (Y CD .GT.0.75  .AND.  Y(I) .LE.C.80)  THEN 
C 

C  Get  the  desired  percentile  at  80% 

C 

DIF80  =  .80  -  Y(I) 

IF  CDIF80.LE.C80)  THEN 
080  =  DIF80 
Y79  =  Y  CD 
DT9  =  D(I) 

Y81  =  YCI+1) 

D81  »  D(I+1) 

ENDIF 

ELSEIF  (Y<I) ,GT, ,80.AND.Y(I) ,LE. .85)  THEN 
C 

C  Get  the  desired  percentile  at  85% 

C 

DIF85  =  .85  -  Y Cl) 

IF  (DIF85.LE.C85)  THEN 
C85  *  DIF85 
Y84  =  YCI) 

D84  =  D(I) 

Y88  =  YCI+1) 

D86  *  DCI+1) 

ENDIF 

5T.SEIF  (Y(I) ,GT, .85.AND.YCI) .LE. .90)  THEN 
C 

0  Get  the  desired  pereen+ile  at  90% 

C 

2IFS0  *  .90  -  YCI) 

IF  CDIF90.LE  090)  THEN 
CPO  -  DIF9G 
Y89  »  YCI) 

D89  =  PCI) 

Y91  =  YCl+i) 

D91  *  DCI+1) 
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ENDIF 

ELSEIF  (Y(I) .GT. .SO.AND.Y(I) .LE. .95)  THEN 
C 

C  Get  the  desired  percentile  at  95% 

C 

DIF95  =  .95  -  Y(J) 

IF  (DIF95  .LE.  C95)  THEN 
C95  =  DIF95 
Y34=  Y  CD 
Y96=  Y(I+1) 

D94=  DCI) 

D96  =  DCI+1) 

EHDIF 

ELSEIF  (Y(I) .GT. .9S.AND.YCI) .LE. .99)  THEM 
C 

C  Get  the  desired  percentile  at  99% 

C 

DIF99  =  .99  -  YCI) 

IF  (DIF99  .LE.C99)  THEM 
C99  ■  DIF99 
Y98  =  Y(I) 

Y100  *  YCI+1) 

D98  *  0(1} 

D100  *  0(1+1} 

EHDIF 

EIDIF 

200  CONTINUE 
300  COMTIMUE 

IF  (DIFSO.EQ.O.O)  THEM 
CRIT80  =  D79 
ELSE 
C 

C  Compute  the  critical  value  at  the  significance  level  of  .20 
C 

CALL  CAH(Y79,Y8l, 079,081,. 80, CRIT80} 

ENDIF 

IF  (DIF8S.EQ.0.0)  THEM 
CRIT85  =  084 
ELSE 
C 

C  Compute  the  critical  value  at  the  significance  level  of  .15 
C 

CALL  CAN(Y84,Y86,D84,D86, .85,CRIT85) 

ENDIF 

IF  (DIF90.EQ.0.O)  THEM 
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CRIT90  =  D89 
ELSE 
C 

C  Coapute  the  critical  value  at  the  significance  level  of  .10 
C 

CALL  CAN(Y89,Y91,D89,D9i, .90,CRIT90) 

EMDIF 

IF  (DIF95.EQ.0.0)  THEM 
CRIT95  =  D94 

‘else 

C 

C  Coapute  the  critical  value  at  the  significance  level  of  .05 
C 

CALL  CAM(Y94fY96»D94,D96, .S5.CRIT95) 

EMDIF 

IF  (DIF99  -EQ.0.0)  THEM 
CRIT99  =  D98 
ELSE 
C 

C  Coapute  the  critical  value  at  the  significance  level  of  .01 
C 

CALL  CAi(Y98,YlOOfD98,DlOO,.99,CRIT99) 

EMDIF 

END 


A-8 


C********************************************************* 


C*  * 
C*  THIS  PROGRAM  EXECUTES  THE  POWER  COMPARISON  * 
C*  FOR  THE  ANDERSON-DARLING  AND  THE  CRAKER-VGK  KISES  * 
C*  TESTS  WITH  TIE  UNIFORM  DISTRIBUTION  AT  THE  GIVEN  * 
C*  SAMPLE  SIZE  * 
C*  * 


O*****************#*******#**********************:****#*** 

n 

V 

C 

c  *****  SOME  COMMENTS  ABOUT  THE  VARIABLES  USED  ***** 

C 

C  N^  Sample  size  used 

C  P0WER(*)  *  The  rejection  number  at  the  given  level 
C  SEED  =  Input/Output  double  precision  variable  assigned 
C  an  integer  value  in  the  exclusive  range 

C  AMBAR,  AMD  *  The  known  and  the  modified  A-D  statistic  values 
C  calculated 

C  CRAKER1 .  CRAMER2  *  The  known  and  the  modified  CV-K  statistic 
C  values  calculated 

C 

INTEGER  N,J,K,L,M,CQUNT(4) ,PQWER(30) *€NTi,€NT2,I 
REAL  WKC360) ,R(120) ,S(120) ,T(120) , 

1  Y , P , ANDAR , AND , CRAMER! , CRAMER2 , PWR ( 30 ) 

DOUBLE  PRECISION  SEED1 
C 

C  Enter  the  double  precision  seed 
C 

SF'D=1095785.D0 

OPat (UNIT=1 ,FILE=’UP80. OUT’ ,STATUS= ’UNKNOWN’) 

DO  tO 0  1=1,30 
?0WER(I)=O 
600  CONTINUE 
C 

C  Enter  the  sample  size 

C 

N=80 

DO  100  .1*1,5000 

y» 

\* 

C  Generate  The  Uniform  deviates 

C 

CALL  GGUBS(SEED1,N.R) 


C 

C 

C 


Sort  the  deviate; 


o  o 


CALL  VSRTA<R,H) 

DO  200  K=1,M 
S(K)  -  R(K) 

200  CONTINUE 
C 

Standardize  the  data 

CALL  ESTPAR(S,N) 

C 

C  Sort  the  standardized  data 

C 

CALL  VSRTA(S,R) 

DO  400  L=1 ,K 
Y=S(L) 

C 

C  Determine  the  CDF  values  of  the  uniform  distribution 
C  for  the  standardized  data 

C 

CALL  fflWORCY.P) 

S(L)=P 
400  COMTIMUE 
C 

C  Find  the  statistic  values  for  the  following  tests 
C 

CALL  AHDARl (M ,  S , ANDAR) 

CALL  ANDAR2(M,S,AMD) 

CALL  CVHKM.S.CRAJERl) 

CALL  CVK2(i,S,CRAMHl2) 

C 

C  Enter  the  critical  values  for  the  sample  size  SO  from 

C  the  critical  value  tables  presented  in  Chapter  4 

C 

IF  (AMDAR.GT. 0.5018578)  THE! 

P0WER(l)sP0WER(l)+l 

EMDIF 

IF  (AMDAR.GT. 0.5535432)  THEM 
P0HER(2)=P0WER(2)+1 
EMDIF 

IF  (AMDAR.GT. 0.6213608)  THEM 
P0WER(3)»P0WER(3)+1 
EMDIF 

IF  (AMDAB.GT. 0.7532310)  THEM 
PQWER(4)=P0WER(4)+1 
EMDIF 

IF  (AMDAR.GT. 1.046329)  THEN 
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P0WER(5)  sPOHER(S)+l 
ENDIF 

IF  CAID.GT.  1. 602745)  THEN 
POHER(6)=PO«LTl(6)+i 
EIDIF 

IF  CAID.GT.  1,650009) THEII 
P0IER(7)-PGWERC7)+1 
ENDIF 

IF  (AND.GT.  I.714il9)TKEl 
PQWER(8)sPQHER(8)+l 
ENDIF 

IFCAMD.GT.  1.824l96)THa 
PQIER(9)SP0WER(9)+1 
EIDIF 

IF(AID.GT.  2. 109204) THEM 
P0HBI(10)=P0WER{10)+1 
EIDIF 

IF CCRAMEE1 , GT . 0 - 0793214445THEI 
?D»3i(ii}=PQiiER(ii)-*i 
SfDIF 

IFCCRAMERl  .GT.0.089231446)THEI 
P0«ER(i2}*POWER(l2)+l 
EIDIF 

IFCCRAMERl .GT.O. 1024222JTHEN 
POIERC 13)*PQHER(  13)+1 
EIDIF 

IFCCRAMERl .GT.O. 12574 16) THEI 
P0WER(14)®P0WER(14)*1 
EIDIF 

IFCCRAIERl  .GT.O.  183783DTHEK 
P0HER(15)=P0HER(15)+1 
EIDIF 

IF (CRAHER2 .GT.O. 083445959)THEI 
POWERC 16) *POWER( 16) + 1 
EIDIF 

IF (CRAMER2 .GT.O. 093562253)THEI 
P0HER(17)SP0WER(17)+1 
EIDIF 

IF  (CRAHER2 .  G7. .  0 . 1065563)7HEI 
P0HERC18) 0WERC18) +1 
EMDIF 

IF (CRAHER2 .67.0.1313581) TKEI 
PGWER(19>=1  n‘47i!.(19)+l 
ENDIF 

IF  (CRAMER2 . ,* . 0 . 1870539)THEH 
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u  o 


PGWER(2Q)=POWER(20}+1 

ENMF 

100  CONTINUE 

DO  500  1=1,20 

Determine  the  rejection  number 
C 

PilR(I}=PQSER(I)/5OO0.O 
500  CONTINUE 

WRITE (1,*) ’FOR  K  =’  ,  K  ,  'POPULATION  =  UNIFORM* 

DO  11  1=1,20 
C 

C  Print  the  rejection  number 

C 

WRITE (1,*)  PHR(I) 

11  CONTINUE 

END 

C 

C  The  subroutines  are  presented  in  the  previous  section 

C 

C 
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